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FARiiUAY  ROTATION 
DIFFERENTIAL  ABSORPTION 
RADIO  WAVE  PROPAGATION 


ELECTRON  DENSITY 
COLLISION  FREQUENCY 


This  report  describes  the  optimization  of  a radio  propagation 
experiment  suitable  for  studies  of  the  ionosphere  D-region  utilizinf 
relatively  low  power,  portable  ground-based  transmitters  and  simple 
receivers  aboard  small  sounding  rockets. 

A wave  propagation  model  has  been  developed  that  numerically 
calculates  the  radio  signal  that  would  be  received  by  a dipole  - 
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'antenna  aboard  a rocket  traveling  up  through  the  ionospheric  D- 
region.  The  Faraday  rotation  and  differential  absorption  experi- 
enced by  the  wave  can  be  used  to  deduce  both  electron  concentra- 
tions and  electron-neutral  collision  rates  in  the  D-region. 

Faraday  rotation  and  differential  absorption  were  numerically 
calculated  for  a number  of  radio  frequencies  for  three  types  of 
extreme  conditions  such  as  occur  at  high  magnetic  latitudes. 

Altitude  profiles  of  these  quantities  provide  Information  to  base 
a selection  of  2 or  3 frequencies  that  will  yield  maximum  infor- 
mation on  electron  density  and  collision  rate,  based  on  the  ex- 
pected ionospheric  conditions.  '>lTie  propagation  experiment  proposed 
employs  multiple  frequencies  tog^her  with  low  power  transmitters, 
and  inexpensive,  movable  antennas Signal  limitations  such  as 
atmospheric  noise,  telemetry  error  and  transmission  power  were  I 

also  used  to  interpret  the  Faraday  rotation  and  differential  I 

absorption  curves  in  order  to  achieve  maximum  accuracy.  .,The  | 

results  of  this  analysis  are:  ] 

1)  For  high  electron  density  (PCA  conditions)  and  low  | 

collision  frequency  conditions,  wave  frequencies  of  I 

3.25,  10,  and  19  MHz  should  bo  used.  j 

2)  For  high  electron  density  and  high  collision  frequency 
conditions,  frequencies  of  3.25,  14,  and  25  MHz  should 
be  used. 

3)  For  low  electron  density  and  low  collision  frequency 

conditions,  frequencies  of  3.25  and  4.25  MHz  are  | 

sufficient.  j 

4)  For  low  electron  density  and  high  collision  frequenev  1 
conditions,  3.25  and  6 MHz  frequencies  are  sufficient.  j 

A method  of  estimating  transmission  power  and  rocket  receiver 
range  requirements  (based  on  noise  levels)  was  developed,  as  well  I 

as  equations  by  which  coning  modulation  of  the  received  signal  can  i 
be  separated  from  actual  Faraday  rotation  and  differential  absorp-  j 
tion.  A data  reduction  scheme  was  developed  to  remove  errors  ' 

associated  witli  the  motions  of  the  rocket.  I 
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CHAPTER  I 
INTRODUCTION 

In  recent  years  several  propagation  experiments  have  been  re- 
ported \_Mcohtly  et  al.,  1967;  Kane,  1959;  Bennett  et  al.,  1972;  and 
tlegill  et  at.,  197l]  which  determine  electron  density  and  electron 
collision  rate  from  the  Faraday  rotation  and  differential  absorption 
of  radio  waves  propagating  between  the  ground  and  a rocket  vehicle. 

It  Is  well  known  that  radio  waves  propagating  In  an  anisotropic  plasma 
medium,  such  as  the  earth's  Ionosphere,  propagate  In  the  Independent 
modes  called  the  ordinary  and  extraordinary  modes  \_Kelao,  1964].  The 
propagation  characteristics  of  each  mode  reveal  that  the  extraordinary 
wave  Is  absorbed  faster  than  the  ordinary  wave,  the  result  being  that 
when  the  propagation  direction  is  near  paralllsm  with  the  geomagnetic 
figure,  the  total  wave,  which  Is  the  sum  of  the  two  modes,  experiences 
an  amplitude  modulation  and  a change  In  the  wave  polarization  as  the 
wave  propagates  through  the  medium.  The  former  effect  Is  called  dif- 
ferential absorption  and  the  latter  Faraday  rotation.  By  measuring 
each  of  these  phenomena  Independently  It  Is  possible  to  work  backwards 
via  a propagation  theory  and  uniquely  determine  the  electron  density 
and  electron  collision  frequency. 

Kane  [1959]  performed  an  experiment  to  measure  electron  collision 
frequencies  when  the  Ionization  density  of  the  ionospheric  D region 
was  anomalously  high.  A 7.75  MHz  radio  signal  and  its  6th  harmonic 
were  transmitted  from  a rocketborne  transmitter  to  a ground  receiver. 

On  the  ground  the  low  frequency  was  multiplied  by  6 and  compared  with 
the  ordinary  component  of  the  high  frequency  signal.  By  utilizing  the 
resulting  beat  frequency  and  the  signal  strengths  of  the  two  polariza- 
tion modes  of  the  7.75  MHz  signal,  it  was  possible  to  determine  the 
collision  frequency  as  a function  of  altitude. 

A high  resolution  propagation  system  to  measure  Faraday  rotation 
and  differential  absorption  was  developed  by  the  University  of  Illinois 
\_Mechtlg  et  at.,  1967].  This  system  employed  an  elaborate  antenna  and 
a constant  power  transmitter  to  launch  a clrcularlly  polarized  3.385  MHz 
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wave  into  the  ionosphere  as  the  ordinary  raode.  The  extraordinary  mode 
was  generated  500  hz  higher  in  frequency,  but  was  part  of  a variable 
power  feedback  loop  which  held  the  extraordinary  signal  strength  at  the 
rocket  (modulation  ratio)  at  a constant  ratio  to  the  ordinary  wave.  By 
carefully  calibrating  the  transmitter  power  and  recording  the  rocket 
roll  rate  by  independent  instruments  it  was  possible  to  extract  both 
differential  absorption  and  Faraday  rotation  data.  With  both  of  these 
propagation  quantities  measured,  both  the  electron  density  and  electron 
collision  frequency  could  be  uniquely  determined  from  the  data. 

f'legill  ct  al . [l97l]  measured  the  electron  density  with  a 3 fre- 
quency (7,  17,  27  MHz)  Faraday  rotation  experiment.  The  data  were  con- 
sidered satisfactory  as  long  as  the  Sen-Wyller  [i960]  and  modified 
Appleton-llartree  \_KeLso,  1964]  forms  (modified  by  using  an  effective 

collision  frequency  \)'  - (5/2)v  where  V is  momentum-transfer  collision 

m m 

frequency)  of  the  wave  propagation  differed  by  less  than  a factor  of  2. 
It  was  found  that  neither  theory  gave  salisfactorv  results  when  the 
radio  wave  frequency  was  less  than  v^/1.2. 

Total  electron  content  along  a propagation  path  was  detennined  by 
Dam  [l972]  by  transmitting  6 continuous  wave  (CW)  frequencies  from  a 
rocketborne  beacon.  Three  frequencies  were  harmonically  related  and 
phase  coherent  so  that  dispersive  phase  measurements  could  bo  made. 

The  other  3 frequencies  were  used  for  determination  of  rf  absorption. 

By  receiving  both  the  left  and  right  hand  polarized  wave  components, 
the  rocket  spin  effeits  and  Faraday  rotation  effects  were  eliminated. 
Electron  density  profiles  could  then  be  determined  from  the  dispersive 
phase'  measurements. 

Hoinett,  7;  . [1972]  reported  a propagation  experiment  where  3 
independent  frequencies  were  transmitted  to  a rocketborne  receiver. 

Each  sign.al  contained  the  effects  of  Faraday  rotation  and  differential 
absorption.  At  any  given  altitude  one  of  the  frequencies  yielded  the 
most  accurate  data,  so  3 values  were  chosen  to  cover  the  altitude  range 
of  interest.  By  utilizing  these  data  with  either  the  Sen-Wyller  or 
Appleton-llartree  magneto- ionic  propagation  theories,  it  was  possible  to 
reconstruct  electron  density  and  electron  collision  rate  profiles. 
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In  view  of  the  complexity  of  some  of  the  above  experiments,  it  was 
desirable  to  investigate  the  possibility  of  a propagation  experiment 
that  would  utilize  simple  and  portable  antennas,  and  inexpensive  trans- 
mitters in  order  that  the  experiment  would  be  modest  in  cost  and  could 
be  transported  to  remote  field  sites.  This  work  describes  such  a mul- 
tiple frequency  propagation  experiment.  The  propagation  model  was 
developed  for  high  magnetic  latitudes  and  includes  the  considerations 
of  errors  due  to  noise,  reflections  and  rocket  coning.  The  signal  re- 
ceived by  a spinning  rocket  was  developed  in  Chapter  II.  These  results 
are  employed  in  Chapter  111  to  illustrate  the  Faraday  rotation  and  dif- 
ferential absorption  that  results  from  extreme  values  of  polar  iono- 
spheric D-region  conditions.  From  these  curves  we  developed  methods 
for  predicting  optimum  wave  frequencies  and  specifications  for  receiver 
gain  levels.  Potential  errors  due  to  rocket  coning  are  evaluated  in 
the  appendix,  and  a method  is  developed  to  systemat i<'al ly  eliminate 
these  errors.  Finally,  a discussion  and  evaluation  of  transmitting 
and  receiving  antennas  is  reported  in  Chapter  IV. 


1 


iT'iCiUi  .a  page  sl..nk-:;ot  p 


s 


CHAPTER  II 
PROPAGATION  THEORY 

The  inital  mathematical  descriptions  of  electromagnetic  wave  pro- 
pagation in  the  ionosphere  were  developed  by  Appleton  [l932]  and  Hartvee 
[l93l].  Their  equation,  which  described  the  refractive  index  of  the. 
wave  as  a function  of  electron  density,  collision  frequency,  wave  fre- 
quency and  propagation  direction  with  respect  to  the  geomagnetic  field, 
has  come  to  be  appropriately  called  the  Appleton-Hartree  equation.  In 
this  derivation  the  collision  process  was  assumed  to  be  independent  of 
electron  energy.  Subsequent  experiments  have  shown  [Phelps  an<f  Pack, 
1959]  that  the  electron-neutral  collision  process  is  indeed  energy- 
dependent.  This  fact  was  utilized  by  Sen  and  Wijller  [i960]  in  an  exten- 
sion of  the  magnetoionic  formulas  of  Appleton  and  Hartree.  Numerous 
experiments  have  shown  this  latter  formulation  to  better  describe  the 
actual  characteristics  of  an  electro-magnetic  wave  in  a plasma  medium. 
Both  of  the  above-mentioned  theories  show  that  a plasma  immersed  in  a 
magnetic  field  is  a birefringent  medium;  i.e.  regardless  of  the  wave 
configuration  impinging  upon  the  medium,  only  two  possible  modes  will 
propagate  internally.  These  propagation  modes  can  then  be  described 

in  part  by  their  complex  indices  of  refraction  n , n and  their  wave 

o x 

polarizations  R^,  R^,  The  subscripts  o and  x refer  to  the  "ordinary" 
and  "extraordinary"  modes,  respectively.  The  refractive  indices  given 
by  Sen  and  WylLer  [i960]  are 

2 / . ,2  A + Bsin^cf)  ± yB^sln**!!  - C^cos^Ji 

n = (m  - JX)  = ^ 

D + Esin^ij) 

where  A,  B,  C,  D and  E are  functions  of  electron  density,  electron- 
neutral  collision  frequency,  electron  gyrof requency , and  wave  frequency, 
and  are  repeated  for  ready  reference  in  Appendix  A.  The  variable  <{>  is 
the  angle  between  the  earth's  magnetic  field  ^ and  the  ray  path.  For 
an  upward  moving  wave  the  quantity  under  the  radical  in  equation  (1) 
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H(z,<}))  = B^sln‘‘4)  - C^cos^(p 


(2) 


is  a complex  number  which  rotates  clockwise  in  the  complex  plane  as 
altitude  z increases.  The  normalized  values  of  equation  (2)  are  repre- 
sented by  the  lengths  and  angles  of  the  arrows  shown  in  Figure  la. 


IMAGINARY 


IMAGINARY 


Figure  la.  Rotation  of  H/|h| 
for  Increasing 
altitude. 


Figure  lb. 


Rotation  of  the 
complex  unit 
TCj  as  altitude 
increases  for 
the  VHF  frequency 
in  Figure  la. 


An  ambiguity  in  the  interpretation  of  the  two  modes  of  (1)  may  arise 
because is  a complex,  double-branched  function  of  a complex  variable. 
The  two  branch  functions  can  be  written 


v/h  = 


Tr,  = 


Tr^  = 


(3) 


where  a Is  the  angle  between  the  real  axis  and  the  radius  vector  de- 
scribed by  the  complex  number  II/|h|.  The  variation  of  Tr^  with 
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altitude  is  shown  in  Figure  lb.  The  ambiguity  in  (1)  can  now  be  removed 
if  Tr,  is  substituted  for  the  radical  in  (1)  giving 


(p  - jx)2  = A + Bsin^(|)  ± Tri 
D + Esin^cfi 


(A) 


The  ordinary  mode  of  wave  propagation  can  now  be  defined  by  using 
the  "+"  sign  in  (A)  and  the  extraordinary  mode  by  the  sign  regard- 
less of  the  relative  magnitudes  of  gyrof requency  and  wave  frequency. 
One  obvious  advantage  of  using  (A)  in  computing  indices  is  that  the 
ordinary  and  extraordinary  components  can  be  determined  independently 
of  the  values  of  gyrof requency , wave  frequency,  collision  frequency, 
and  plasma  frequency,  which  are  imbedded  in  the  angle  a.  In  addition, 
the  discontinuity  encountered  by  Ratcliffe  [l962](pp.  69-70  and 
Figure  7. A)  no  longer  occurs.  The  function  H also  appears  in  the  wave 
polarization  equation  of  Sen  and  Wyllev  [i960]  which  can  be  written  in 
the  form 


Bsin^c))  ± Trj 
Ccos(^ 


(5) 


where  again  + and  - correspond  to  the  ordinary  and  extraordinary  wave 
modes,  respectively. 

A coordinate  system  (jt,  is  defined  with  z vertical  and  2. 

pointing  magnetic  north.  A typical  rocket  trajectory  illustrated  in 
Figure  2 is  described  in  this  system.  Only  the  trajectory  upleg  will 
be  considered  in  this  report.  The  wave  propagation  is  defined  in  the 
(ji'.  X.'  > coordinate  system,  where  z'  is  parallel  to  the  direction 
of  propagation  and  lies  in  the  magnetic  meridian.  The  angle  <)  is 
the  angle  between  the  geomagnetic  field  in  the  northern  hemisphere  and 
the  upward  ray  path.  The  wave  polarization  in  the  latter  coordinate 
system  is  given  by 


R 

o 


ox 


+ 31, 


with  <S  ,Y  real 
o o 


o 


(ha) 
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-E 

XX 


where  the  total  ordinary  wave  component  of  the  electric  field  is  E = 

— o 

E x'  + E y'  and  the  extraordinary  component  is  E = E x'  + E y' 

ox—  Oy^  J r ^ 

The  polarization  is  defined  by  the  ratio  of  the  transverse  electric 

field  components.  In  general,  both  modes  are  elliptical ly  polarized 

since  a complex  value  of  R means  elliptic  polarization  [Kelso,  1964, 

pp.  36,37].  However,  for  our  geometry  and  for  frequencies  typically 

of  interest  in  Faraday  rotation  and  differential  absorption  measurements, 

i.e.  between  2.5  and  25  MHz,  the  values  of  6^  and  ) ^ are  approximately 

10“^  while  0.998  < y 1.00  and  -1.00  < v < -0.998.  Thus,  for  the 
o X 

geometry  of  Figure  2,  the  wave  polarizations  of  equation  (6),  can  be 
approximated  to  good  accuracy  by  R^  = +j  and  R^  = -j  on  the  upleg  of  the 
rocket  flight.  R = ij  corresponds  to  circularly  polarized  waves  in  the 
quasi-longitudinal  propagation  approximation.  Figure  3 shows  schemati- 
cally that  R = +j  and  R = -j  . This  means  that  in  the  Northern  hemi- 
sphere,  the  ordinary  mode  ^-field  rotates  clockwise  in  time  when  viewed 
along  the  +z'  direction,  while  the  extraordinary  mode  vector  rotates 
in  the  opposite  direction.  Actually,  the  more  general  elliptical  waves 
of  (6)  also  rotate  in  the  clockwise  and  counterclockwise  directions. 

We  assume  that  the  initial  wave  is  excited  with  the  electric  field 
aligned  parallel  to  the  x'  axis  (Figure  3).  Below  the  D-region,  this 
^-field  can  be  described  by 


E 


Mej(‘^t  - kz')^. 


(7) 


where  k is  the  propagation  coefficient,  to  is  the  wave  angular  frequency 
and  M is  the  wave  amplitude.  For  R ~ ±j , the  wave  described  by  (7)  can 
be  written  as  the  sum  of  two  circularly  polarized  waves,  an  ordinary 
mode 


M j(u)t  - k z')  , .M  j(u)t  - k z')  , 
o X ■ 3"^  o y 


2® 


E 


(8a) 
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Figure  3.  EM  wave  coordinate  system.  ^ lies  in  the  , ^')  plane. 

Propagation  is  in  the  (^')  direction  with  rotation  of  the 
two  modes  as  indicated. 


and  an  extraordinary  mode 


E 

— X 


M j (cjt  - k z ' ) , 
X x' 


+ 12® 


(cot  - k z')  , 
x 1 


(8b) 


Since  the  Faraday  rotation  and  differential  absorption  will  be  shown  to 
depend  only  on  the  difference  between  the  two  allowed  modes  of  propaga- 
tion, the  wave  phase  at  the  bottom  of  the  ionosphere  is  arbitrary;  for 
con’-enlence  we  assume  it  to  be  zero.  Note  that  (8a)  and  (8b)  add  to 
give  the  original  plane  wave  (7).  The  two  propagation  coefficients  are 
given  by 


k 

o 


w w,  , . 

- n = -(u  - jY  ) 

CO  C O ^o 


(9a) 
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, w w,  . . 

k = — n = — (u  - 1 Y ) 
X CO  c X ''■X 


where  quantities  Xq>  depend  on  electron  denstiy  N^,  electron 
collision  frequency  V^,  the  earth's  magnjitic  field  ^g,  and  the  angle  of 
propagation  (]).  If  and  vary  with  altitude  only,  then 


Uq  = Vig(z,<j>) 


Xo  = Xg(z,<})) 


Mx  = U^(z.(D) 


(10b) 


Let  0(z)  be  the  angle  between  the  ^ and  z'  axes  of  Figure  2.  The  path 
length  is  then  given  by  z'  = z sec[0(z)]  where  z is  the  vertical  height. 
By  substituting  (9)  and  (10)  into  (8) , the  electric  field  components  of 
the  two  ionospheric  modes  at  a distance  h'  into  the  plasma  can  be  written 


^(h')  = |exp  j j[a)t  - 4)^(h')]  - ^^(h’)  j x' 


+ jfexp  ' j[u)t  - (t»  (h')]  - ot  (h’)  [ x' 


E (h')  = ^xp  j[tot  - <{)  (h')]  - a (h')  x' 

X Z ( O X I 

- jyexp  ) j[ajt  - (j>^(h’)]  - a (h')  [ 


where 


•'o 

!>x(h')  = f / M^(z,0)dz’ 


- / u^(z,4))dz’ 


(12b) 
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ao(h') 


a^(h') 


h' 

X^(z,(t))dz' 


h' 

X^(z,<J))dz' 


(12c) 


(12d' 


It  is  shown  in  Appendix  B that  the  curvature  of  ray  paths  in  the  D-region 
is  extremely  small;  therefore,  it  is  assumed  that  0(z)  is  the  same  for 
both  the  ordinary  and  extraordinary  waves  and  that  6 remains  constant 
over  a given  path.  At  altitude  h,  h = h'cos[0(z)],  the  total  E-field  is 
given  by  the  vector  sum  of  (11a)  and  (lib).  The  observable  part  of  the 
total  electric  field  would  then  be 


I -a  (h) 

Re[j%(h)]  ~ J I ^ cos[wt  - ({)^(h)]  + e cos[oot  - ())^(h)] 

I -a  (h)  -a  (h) 

- j e ° sin[wt  - ^’^(h)]  - e ^ sin[uJt  - <!)^(h)] 


Z'  . (13) 


At  a given  altitude  h.  Re  ^(h)  represents  a vector  with  direction- 

dependent  magnitude  rotating  in  space  with  frequency  O).  The  locus  of 

points  traced  out  by  (13)  is  an  ellipse.  The  maximum  length  of  Re(E) 

will  not  necessarily  be  parallel  with  either  the  >c'  or  axes,  nor 

are  the  values  of  E ' and  E ' the  same. 

X y 

The  next  step  in  the  interpretation  of  the  wave  propagation  is  to 
transform  Re[E(h)]  of  (13)  to  a coordinate  frame  which  rotates  in  syn- 
chronism with  the  polarization  ellipse.  Since  4)  and  d)  are  functions 

o ’^x 

of  heiglit,  the  orientation  of  the  ellipse  rotates  with  altitude.  The 
transformation  from  the  electromagnetic  wave  coordinate  axes  (x ' , y' , 
^')  of  figure  3 to  the  axes  (x*^,  ^^)  rotating  about  the  axis  is 

given  by  the  matrix 


R(h) 


cos  [(4>^  - <P^)/2  ] 
sin  [(4)  - 4)  ) /2  1 

O X 


sin  [(4)^  - <P^)/2  ] 

cos  [(4)  - 4)  )/2l 

o X 


0 

0 


. (14) 


0 


0 


1 
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R is  seen  to  be  orthogonal  [Det(R)  = l],  so  that  R represents  a proper 
rotation  from  the  (x' , ) axes  about  the  z'  axis  of  ((|)^  - (})^)/2  radians. 

If  (13)  is  written  as  a column  vector 


Re(E)  = 2 


-a  -a 

e cos  (cut  - <p  ) + e cos  (cut  - <]>  ) 
o 

-u  -a 

-e  sin((ut  - i|)^)  + e sin(iut  - <))^) 


0 


(15) 


the  matrix  product  R • [Re(E)]  yields  the  new  vector 


'~\X  ~(X 

(e  + e ^)cos|tiut  - ((})^  + 4)^)]/2| 

^)sin|LijJt  - (4)^  + 4>^)]/2| 


-u  -a 
(e  “ - e 


(16) 


One  can  easily  verify  that  in  the  (x*^,  reference  frame 

represents  the  equation  of  an  ellipse  with  major  axis  parallel  to  x*^ 

and  minor  axis  along  The  relationship  is  illustrated  schematically 

in  Figure  4 where  the  angle  (<{>^  ~ angle  of  rotation  between 

the  (x’ , y^')  and  (jc"^ , y*^)  coordinate  axes.  The  lengths  of  the  major  and 

M ~oi  -a  M -a  -a 

minor  axes  are  respectively  + e ^)  and  ° - e ^) . The 

rotation  between  the  two  coordinate  systems  (x' , y')  and  (x’^,  y*^)  is 
called  Faraday  rotation  with  F = ((})^  - 4'j^)/2  being  the  Faraday  rotation 
angle  and  the  shape  of  the  ^-field  ellipse  vary  with  position  along  the 
rocket  trajectory.  Calculations  show  that  i})^  > (J)^  for  all  measurable 
quasi-longitudinal  cases,  so  not  only  does  rotate  clockwise  within 
the  ellipse  at  frequency  (u,  but  the  ellipse  itself  rotates  slowly  clock- 
wise as  altitude  increases,  as  described  by  the  Faraday  rotation. 


Figure  4.  Relation  of  the  x’",  ^ axes  to  the  x',  axes  for 
two  different  altitudes,  where  h,  ^ hj. 


Differential  absorption  A is  defined  as 


A = 20  log 
^ 1 0 


( 


In  tile  fr.inie  the  electric  field  amplitude  can  be  written 


E 

o 


E,  + 


for  the  ordinary  wave,  and  for  the  extraordinary  mode 
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where 

(19a) 

(19b) 

Substitution  of  (18)  and  (19)  into  (17)  yields  an  expression  for  the 
differential  absorption. 

A = 20[a^  - a^llogj^e  (20) 


As  discussed  earlier  both  the  Faraday  rotation  and  differential 
absorption  (14)  and  (20)  are  independent  of  the  initial  phase  and  ampli- 
tude of  the  wave  at  the  bottom  of  the  ionosphere. 

It  is  now  assumed  that  the  wave  transmitter  is  physically  near  the 
rocket  launch  site,  and  the  rocket  symmetry  axis  ^ which  is  parallel 
to  the  direction  in  which  the  rocket  is  initially  fired,  is  constant 
over  the  up-leg  of  the  rocket  flight.  Then  and  are  almost  par- 
allel, see  Figure  2.  For  a rocket  spinning  in  the  same  sense  as  the 
Faraday  rotation,  the  transformation  from  the  frame  in  which  the  rocket 
antenna  is  stationary  (jt^,  to  the  electromagnetic  wave  frame 

(2<' , x' ’ Z.' ) represented  by 


S = 


cos  at 
sin  at 
0 


-sin  at  0 

cos  at  0 


0 


1 


(21) 


where  a is  the  spin  rate.  The  transformation  from  the  (x'^ , , ^‘*)  axes 

to  the  Faraday  rotation  axes  is  given  by  (14)  and  (21) 
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T = RS  = 


cos [at 

-(({,- 

4>  )/2] 

-sin[at  - 

(<t> 

- 

)/2] 

o 

X 

o 

X 

sin[at 

-(>!>- 

4>  )/2] 

cos[at  - 

((p 

- p 

)/2] 

o 

X 

o 

X 

0 

0 

rocket 

antenna 

to  be  a 

short  dipole 

pointing 

rad 

(22) 


ward,  i.e.  perpendicular  to  the  rocket  axis.  For  other  antenna  configur- 
ations, this  would  be  the  direction  of  the  equivalent  dipole.  Let 


T 


a 


1 

0 

0 


represent  a unit  vector  in  the  dipole  direction,  i.e.  x^  is  parallel 
to  the  antenna.  Multiplying  on  the  left  by  T gives 


= Tt 


a _ 


cos[at  - ((p^  - <p^)/2] 

sin[at  - (i}*  - (|)  )/2] 
o X 

0 


(23) 


which  represents  the  rocket  antenna  direction  in  the  Faraday  rotation 
frame  (x^,  y'',  z’’).  The  relationship  between  the  various  coordinate 
axes  is  shown  In  Figure  5.  Now  that  the  transformations  between  the 
three  principal  coordinates  are  available,  the  value  of  the  ]l-field  in 
the  direction  of  the  rocketborne  antenna  can  be  found  by  forming  the 
scalar  product  of  (16)  and  (23): 


E 

s 


M 

2(e 


-ex 

o ) cos  [cot 


(({)  + <p  )/2]cos[at 

o x 


(<-0  - 


M 

2 


-a 

^)sin[ait 


(ip  + (p  )/2jsin[at 
o X 


($^^  - <,\)/2]  (24) 


e 


/ 


I 

y 


Ki>;uro  5.  Ro  1 ;U  ionshi p botwi’on  coordinate  axes. 


Tile  potential  measnreii  at  the  antenn.i  terminals  is  sim|ily 


wliere  t is  the  effective  electrical  length  of  the  equivalent  aipole 

antenn.i.  Sp.it  ial  variations  transverse  to  the  directions  of  propagation 

can  be  neglected  if  we  use  simple  horizontal  dipole  antennae  at  the 

transmitter.  From  (24)  it  can  be  noted  that  V is  a scalar  function  of 

s 

time  t,  spin  rate  (t,  and  altitude  z through  the  changes  in  a^,  Oi^, 
and  with  electron  density  and  collision  rate.  When  tlie  altitude  is 
held  fixed,  tlu-re  is  a r.ipid  oscillation  in  at  frequencyd'  with,  a much 
slower  amplitude  modul.it  ion  a.  This  is  exactly  what  would  be  expected 
from  a rot.itirig  dipole  antenna  in  the  field  of  .i  plane  electromagni't  ic 
w.ive . 

For  a F.irad.iv  rot.it  ion-d  i f f erent  ia  1 absorption  measurement,  the 
values  of  uj  = 10  sec  .and  a ~ 20  tt  sec 


might  be  typiial.  This  means 
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that  the  rocket  antenna  does  not  move  significantly  during  the  period 
of  one  wave  cycle.  The  voltage  waveform  at  the  antenna  terminals  then 
consists  of  a sinusoid  signal  of  angular  frequency  lO  ~ 10^  sec  with 
a much  slower  amplitude  variation  produced  by  rocket  spin  and  Faraday 
rotation  not  exceeding  10^  sec  The  maximum  value  of  the  induced 

signal  voltage  V during  one  cycle  of  the  wave  frequency  (l  = Zir/to 
= b X 10”^  sec)  is  proportional  to  the  value  of  when  dE^/dt  = 0. 
Differentiating  (24)  with  respect  to  time,  but  holding  at  constant  be- 
cause the  spin  is  quasi-stationary  with  respect  to  o),  gives 

= -Eja)sin[wt  - (cj)^  + cf)^) /2]cos[at  - (({)^  - ({)^)/2] 

-Ej,a)cos[a)t  - (ip^  + (f)^) /2]sin[at  - - (p^)/2]  . (26) 

Setting  (26)  equal  to  zero  gives  avalueof  tot  - (A  + 6 )/2  for  which 

o X 

E is  maximum,  i.e. 
s 

[wt  - (iP^  + = -tan"'  I tan[at  - (4)^^  + 4>j^)/2]|  (27) 


Ej  = M/2(e  ” + e ^) 

-a  -a 

= M/2(e  “ - e ''')  . 

Substituting  (27)  into  (25)  yields  the  peak  value  of  V over  a wave 

s 

period 


''r  = ^ 

1 -1  1 
cos (tan  ; — 

\ < ‘•'1 

tan[at  - 

|cos[at  - 

(^o 

- 

+ V 

2 

sin(tan-'^  -- 

tan[at  - 

(^o 

1 sin[at  - 

(^o 

- 

>■ 
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where  = R.E^,  • Passing  the  instantaneous  antenna  voltage  (25) 

through  a detector  with  a time  constant  long  compared  with  the  wave 
period,  the  detected  output  voltage  would  be  given  by  (28).  This  vol- 
tage contains  the  differential  absorption  and  Faraday  rotation  informa- 
tion about  the  ionosphere  and  needs  only  be  interpreted  to  evaluate  the 
D-region  electron  density  and  collision  frequency. 

The  differential  absorption  given  by  equations  (1/)  through  (20) 
requires  values  for  E^/E^.  These  are  proportional  to  the  maximum  and 
minimum  values  of  V which  occur  when  the  antenna  is  parallel  to  the 
major  and  minor  axis  of  the  ^-fieid  ellipse.  Figure  6a  illustrates  the 
E-field  polarization  ellipse  geometry.  Since  V is  proportional  to  the 
maximum  value  of  (^^  • ^’’)  over  each  cycle,  it  is  seen  that  when  the 
major  axis  of  the  polarization  ellipse  is  parallel  or  perpendicular  to 
the  antenna  the  value  of  V is  maximum  or  minimum.  Another  aspect  of 
the  signal  is  illustrated  in  Figure  6b  where  the  time  varying  voltage 

is  displayed.  The  envelope  of  this  signal  is  V and  the  minimum  and 

R 

maximum  values  of  V_  are  V,  and  V..  In  reality,  the  oscillation  of  V 

R 1 z ^ ’ s 

is  much  more  rapid  relative  to  oscillation  of  V . Since  the  shape  of 

R 

the  polarization  ellipse  changes  as  the  wave  propagates  through  the 
ionosphere,  the  values  and  Vj  also  change  and  the  determination  of 
differential  absorption  depends  upon  their  evaluation.  Faraday  rota- 
tion is  also  obtained  from  the  voltage  of  Figure  6b,  but  it  is  more 
easily  understood  after  several  graphs  are  developed  in  Chapter  III. 


Slil.KCTION  OK  KRKQUKNCIKS 


In  this  chapter  we  use  the  theory  ot  Chapter  11  to  develop  a scheme 
lor  selecting  wave  Ireqnencies  that  will  yield  optimum  results  t roin  a 
I)-region  Karaday  rotation  - d i f lerent  ia 1 absorpl ion  experiment.  A com- 
puter program  was  developed  to  calculate  numerically  the  radio  wave 
signal  that  would  be  detected  at  a rocket  vehicle  traveling,  through  the 
ionosphere  on  a given  trajectory.  A model  ionosphere  ol  the  electron 
density  and  e lec  t ron-neu  t r.i  1 collision  rate  is  required  ti'  calculate 
wave  prop.iga t ion . Profiles  used  for  illustration  were  chosen  to  bracket 
extreme  conditions  that  might  be  found  in  the  polar  l)-region.  Electron 
density  profiles  for  "quiet"  and  "disturbed"  conditions  were  measured 
during  PCA  events  and  are  shown  in  Figure  7.  These  conditions  are  not 
unique  to,  but  may  be  typical  of  rather  severe  PCA.  Collision  frequency 
versus  altitude  is  plotted  in  Figure  8.  The  winter  profile  of  collision 
rate  is  a moderate  estimate  that  might  be  expected  over  the  winter  months. 
At  high  l.ititudes,  the  summer  season  and  possibly  other  factors  can  lead 
to  large  variations  in  the  collision  frequency  [)rofile,  as  commented 
upon  by  .’iii'anr  et  a.'.,[l96b]  and  BclrotiP  el  a/.,[l9bb].  We  have  chosen 
a high  collision  rate  of  three  times  the  lower  curve  to  bracket  the 
possible  values. 

The  profiles  of  disturbed  D-rcgion  electron  density  and  winter 
collision  rate  were  used  to  calculate  a rocket  antenna  signal.  Figure  9 
illustrates  a 0.5  second  segment  of  (equation  28)  each  5 km  altitude 
for  a wave  fre(|uency  of  7 MHz.  The  wave  at  the  bottom  of  the  ionospliere 
was  assumed  to  be  linearly  polarized,  so  that  the  minimum  value  of  V 

R 

(V^  in  Figure  6b)  is  essentially  zero  at  50  km.  As  the  wave  propagates 
higher  into  the  ionosphere,  the  polarization  ellipse  gets  "fatter"  with 
the  result  that  the  minimum  value  of  V gets  larger,  as  illustrated  in 
Figure  9b  and  9c.  Similar  curves  for  20  MHz  are  plotted  in  Figure  10. 

The  liigher  frequency,  20  MHz,  is  less  effected  by  tlie  ionosphere  than 
the  7 MHz  signal.  Since  differential  absorption  depends  upon  the  largest 
and  smallest  values  of  Figures  9 and  10  illustrate  liow  differential 


r 


ELECTRON  DENSITY  (cm“^ ) 


Kiv’iiro  7.  i;  I ec  t ron  donsity  profiles.  Tlie  profiles  labeled 

* QU 1 bl  * and  * D1  STIIRRKI) ' wore  chosen  by  the  authors 
to  bracket  a wide  range  of  possible  conditions. 


L. 


ALTITUDE  (km) 


disturbed  electron  density  and  winter  collision  frequency 
et  altidues  of  60  km  (top)  and  65  km  (bottom). 
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profiles  for  rocket  altitudes  of  70  kni  (top)  and  75  km  (bottom).  Note  the 
potential  scales  have  been  multiplied  by  factors  of  0.1  and  0.01. 
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Figure  lOci.  20  MHz  signal  for  disturbed  electron  density  and  winter  collision  frequency 
profiles  for  rocket  altitudes  of  85  km  (top)  and  95  km  (bottom). 
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absorption  may  bo  extracted  from  tlie  telemetry  data,  'i'be  lowest  altitude 

at  which  differential  absorption  can  be  measured  is  the  lowest  altitude 

where  V can  be  identified  to  be  at  its  minimum  values.  This  measure- 
K 

ment  is  limited  by  noise,  as  discussed  latter. 

The  voltage  V given  by  equation  (28)  is  composed  of  two  terms;  each 

K ^ 

term  is  a periodic  function  of  the  "spin"  argument  [cit  - (4)^  - ({)^)/2]  . 

The  maximum  values  of  V„  occur  at  [at  - (<p  - <p  ) /2]  = nil  and  minimum 

R o x 

values  at  [at  - (i})^  - (})^)/2]  = (n  + 1/2)tt,  where  n is  an  integer.  The 

Faraday  rotation  [(0  - <p  )/2]  can  be  deduced  from  the  periodicity  of  V , 

^ o X R 

if  a is  known.  If  (({>  - 4>  )/2  were  constant  over  a number  of  cvcles  of 

o x _ ' 

at,  then  V would  obtain  its  minimum  value  V every  t = n/a  seconds  cor- 
responding  to  a half  rotation  of  the  rocket.  On  the  other  hand,  if 

varies  slightly  during  each  period,  then  the  minimum  value 
would  be  slightly  displaced  in  time  from  a regular  period.  A rocket 
spin  rate  of  10  rps  (20/71  rad/sec)  was  used  in  calculating  the  data  of 
Figures  9 and  10.  Tick  marks  every  half  spin  revolution  are  marked  on 
the  lower  boundary.  These  marks  coincide  with  the  antenna  lying  parallel 
to  the  magnetic  meridian  plane.  It  can  be  seen  that  the  signal  minima 
and  the  tick  marks  shift  apart  as  the  wave  propagates  upward.  A measure 
of  this  Separation  determines  the  Faraday  rotation.  In  an  actual  experi- 
ment, the  rocket  spin  rate  must  be  resolved  by  an  independent  measurement 
such  as  an  on-board  magnetometer.  The  spin  rate  data  can  then  be  used  to 
help  recover  the  Faraday  rotation  rate. 

Radio  interference  at  frequencies  of  interest  for  a rocket  experi- 
ment comes  fia'm  other  radio  frequency  transmissions,  from  atmospheric 
and  cosmic  noise,  and  from  receiver  and  telemetry  equipment.  Figure  11 
shows  the  relative  magnitude  of  some  of  these  noise  sources.  The  root 
moan  square  electric  field  corresponding  to  a given  noise  power  is  pro- 
portional to  an  equivalent  antenna  noise  temperature  T^  and  bandwidth  B 
by  the  relationship 

= [l20nBkT^]^  (volts/meter) 


where  k is  Boltzmann's  Constant  (1.38x10  ^^watt-sec/K  ).  On  the  basis  of 
(29)  we  can  talk  about  the  equivalence  of  noise  temperature  and  noise  power. 


Antenna 


Fi);iiro  11.  Noise  torapcrnturc'  .it  mocliiim  .inci  hip.h  t requpne  ies 
[/  .F  , 1968].  Curve  Inbeleil  'Test 

noise'  w;is  selected  hy  authors  for  oil  rompvilnt ions 
involving 
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One  other  important  factor  in  receiver  design  is  the  exclusion  of 
unwanted  radiofrequency  transmissions  [Ointher  and  Smith,  1975].  A 
foreknowledge  of  transmissions  at  the  experiment  site  allows  one  to  se- 
lect frequencies  close  to  those  listed  above,  but  with  a minimum  of  inter- 
ference. Such  a choice  together  with  a narrow  receiver  bandwidth  minimizes 
the  error  from  this  source.  Some  serious  consideration  should  be  given 
to  this  possible  source  of  noise  pollution  before  the  experiment  para- 
meters are  firmly  fixed. 

The  total  noise  can  be  represented  as  an  equivalent  electric  field 
and  can  be  written  as  the  sum  of  noise  generated  outside  the  rocket,  sig- 
nal degradation  by  the  telemetry  system  and  receiver  generated  noise 


^REC 


(30) 


In  the  1-30  MHz  region,  a good  receiver  has  a noise  figure  of  about  3 dB, 

which  means  the  receiver  equivalent  noise  temperature  is  only  twice  that 

of  an  ideal  noiseless  receiver.  This  means  that  the  receiver  noise  can 

generally  be  neglected  compared  to  the  other  noise  sources.  Telemetry 

noise  is  expressed  as  a small  fraction  S of  the  full-scale  signal  E,,  ; 

K K j 

the  full-scale  voltage  recorded  on  the  lowest  gain  range  of  the  receiver. 

A good  telemetry  system  will  have  a value  of  S around  0.01.  The  receiver 

K 

system  should  be  designed  so  that  the  telemetry  error  is  approximately 

the  same  order  of  magnitude  as  the  other  noise  sources.  If  >>  E,, 

R R j A 

then  the  telemetry  degradation  is  unnecessarily  large  while  if  E.  >>  S E 

A R R j 

the  receiver  may  be  unnecessarily  expensive.  Assuming  that 


Sr^r,=  ^-A 


(31) 


the  noise  equivalent  electric  field  is  approximately 


E = 2E, 
N A 


If  we  ignore  signals  that  are  less  than  10  dB  above  the  noise  level. 
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the  minimum  signal  that  would  yield  accurate  data  is 

E,  = 2v/l0  E^  . 


(32) 


In  terms  of  ground-based  transmitter  power  P^,  transmitter  antenna 
gain  G^,  and  the  absorption  Indices,  (32)  becomes 

= 2/10  E R (30  P G )■'"  (33) 

Am  11 

where  h^  is  the  lowest  altitude  at  which  condition  (32)  is  satisfied  and 

R is  the  distance  from  the  transmitter  to  the  rocket  receiver.  This 
m 

minimum  altitude  will  obviously  depend  on  transmitted  power.  Thus,  the 
stronger  the  transmitted  signal  the  lower  the  altitude  and  electron 
density  that  can  be  measured  by  differential  absorption. 

The  maximum  received  signal  is  proportional  to  Ej , equation  (19). 

At  the  altitude  where  (33)  is  valid  E^  can  be  written 


-a„(h  ) -a  (h  ) 
o m X m 

e - e 


Ej (h^)  = (30  P^  G^) ^ 


-a  (h  ) -a  (h  ) 
o m X m 

c + e 


R ■* 
m 


(34) 


The  difference  between  E and  E,  at  h is  several  orders  of  magnitude. 

1 2 m 

Therefore,  it  is  necessary  to  incorporate  an  automatic  gain  switching 
capability  into  the  rocketborne  receiver  in  order  to  measure  linearly 
this  large  variation  in  signal  strength.  This  means  that  several  ranges, 
each  liaving  a different  value  of  gain,  must  be  used.  If  each  range  has 
an  impedanct!  separating  it  by  15  dU  from  the  next  highest  range,  tlien 
the  k^^  range  can  measure  a field  strength  of 


E 


RK 


The  range  that  measures  the  largest  voltages  has  the  lowest  gain  and  must 

measure  E (h  ) linearly.  For  S = 0.01,  E = E /O.Ol  and  this  implies 
* m K K ^ A 

that  the  largest  value  of  k is  the  integer  < given  by 
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4 

K > ^ log 


1 0 


E,  (h  ) 
1 m 


- 1.667 


(35) 


Thus  K is  the  number  of  ranges  separated  by  15  dB  necessary  to  measure 

E, (h  ) linearly. 

' m 

The  necessary  criteria  to  estimate  the  lowest  altitude  at  which 
differential  absorption  measurements  can  be  made  are  now  available. 

Some  compromise  will  probably  be  necessary  in  the  choice  of  transmitted 
power,  number  of  receiver  ranges,  etc.,  but  there  is  enough  flexibility 
to  satisfy  most  investigations. 

There  also  exists  a practical  upper  altitude  limit  for  measurement 
of  differential  absorption.  This  limit  arises  from  the  progressively 
decreasing  separation  between  signal  maxima  and  signal  minima,  Ej  and  E^ . 
Figure  9c  illustrates  the  antenna  voltage  at  an  altitude  where  the  extra- 
ordinary wave  has  been  almost  completely  absorbed.  At  some  altitude  h^, 
Ej  and  E^  will  become  indistinguishable  due  to  noise  fluctuations.  For 
the  D-region  the  ordinary  wave  component  is  usually  so  large  that  the 
signal  must  be  measured  in  one  of  the  middle  ranges  of  the  receiver. 

Here  the  noise  is  caused  primarily  by  telemetry  degradation.  Thus 


E..  = S,.  E 


R Rk 


= S. 


10  ^ E„ 


If  we  require  Ej  - E^  to  be  10  dB  above  the  noise,  then 

E,(h^)  - E^(h^)  > 10\ 

Since  Ej^j^  > Ej(h^)  and  lO^'^'^Ej  (h^)  > Ej^j^  = Ej^/S^^,  the  highest  altitude 
at  which  differential  absorption  may  conservatively  be  evaluated  is 


^(h  ) - 


E,(h^)  > 


lo'^Sj^  10 


3/4 


10"s^ 


E /S 
N R 


By  substituting  E^ 


E - E and  E, 
o X 1 


E + E in  the  above  equation,  the 
o X 
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following  expression  results 


2E  (h  ) > 10^-^^S„[E  (h  ) + E (h  )] 

XZ  Roz  K Z ■“ 


whicli  can  he  rearranged  to  give 

(h  ) 


-20.2  dB 


(36) 


if  = 0.01.  Equation  (36)  provides  a means  to  determine  the  approximate 
altitude  where  the  differential  absorption  becomes  too  small  for  measure- 
ment. Equation  (20)  was  integrated  nvin^er ically  for  the  conditions  of 
Figure  9.  The  resulting  differential  absorption  is  shown  in  Figure  12. 

The  lowest  altitude  at  which  differential  absorption  can  be  measured,  the 

condition  of  (33),  was  computed  for  several  different  values  of  T . De- 

A 

creasiii!’  transmitter  power  or  increasing  Llie  bandwidth  by  an  order  of 

magnitude  will  have  the  same  effect  on  the  lower  limit  for  differential 

absorption  ;is  increasing  the  antenna  temperature  by  an  order  of  magnitude. 

Based  on  the  recommendation  of  Westlund  (Chapter  IV),  values  of  = 

200  watts  and  = 3.5  were  chosen  for  this  example.  A bandwidth  of 

2 Khz  was  chosen  for  this  illustration.  For  a median  value  of  noise 

T^^  = 2 X 10‘‘K^,  See  Figure  11,  E^  = 0.9  7 x lO^^’V/m.  The  minimum  height 

for  differential  absorption  measurements  is  calculated  to  be  h - 48.1  km: 

111 

tile  maximum  electric  field  at  this  altitude  is  Ej(hj^)  ~ 6.5  x 10  ^V/m. 
Therefore,  tlie  minimum  and  maximum  antenna  voltages,  and  V^,  differ 

"R 


by  three  orders  of  magnitude.  From  (31)  E ~ 1 x 10  ‘*V/m,  and  condi- 


tion (35)  specifies  4 receiver  ranges.  The  full  scale  equivalent  elec- 
tric field  strength  for  each  receiver  range  and  the  equivalent  telemetry 
noise  are  given  in  Table  I. 

The  -20.2  dB  limit  of  (36)  occurs  at  72.9  km.  Figure  12.  Above  this 

altitude  the  rate  of  change  of  the  differential  absorption  is  masked  by 

the  telemetry  uncertainty.  At  72.9  km,  E (ii  ) - E (h  ) - 0.4  x 10”^V/m 

I z ? z 

which  would  be  measured  in  the  third  range.  Table  I indicates  that  the 
telemetry  error  is  a factor  of  10  larger  than  atmospheric  noise  in  this 
receiver  range. 


ALTITUDE  (km) 
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Curves  similar  to  that  of  Figure  12  can  be  calculated  for  a variety 
of  frequencies.  The  limitations  on  differential  absorption  measurements 
can  be  imposed  and  the  set  of  curves  used  to  predict  appropriate  fre- 
quencies which  will  give  good  results  and  cover  the  desired  altitudes. 
These  curves  are  based  of  course  on  certain  expected  electron  density 
and  collision  rate  profiles. 

TABLE  I 


Err (mV /m) 


s,  E„^(pv/.) 


The  limitations  associated  with  determination  of  electron  density 
by  measurement  of  Faraday  rotation  are  partly  due  to  the  differential 
absorption  phenomenon  and  partly  due  to  telemetry  bandwidth.  The 
highest  altitude  of  measurement  is  reached  when  the  extraordinary  wave 
is  completely  absorbed.  This  limit  is  then  the  same  as  the  differen- 
tial absorption  limit.  The  lowest  altitude  at  which  Faraday  rotation 
can  be  measured  is  determined  by  the  ability  to  detect  a shift  in  the 
signal  null  position  from  the  spin  null  position.  Limiting  the  receiver 
bandwidth  also  limits  the  ability  to  distinguish  the  precise  time  of  the 
signal  nulls.  A 1°  resolution  was  achieved  by  Knoehel  arid  Skaperdas 
[l966].  This  value  will  be  adopted  for  our  determination  of  Faraday 
rotation.  We  note  th.it  the  width  of  the  Faraday  rotation  minimum  in- 
creases as  the  separation  between  signal  maxima  and  minima  decreases, 
i.e.  with  incre.asing  altitude.  At  the  lowest  altitudes  •'‘he  minima  are 
reason.ibly  sharp  and  will  permit  Faraday  rotation  measurements  to  be 
sensitive,  even  though  the  rate  of  change  of  the  phase  is  not  very  largo. 

Actually,  the  m.agnetoionlc  theory  of  wave  propagation  is  very  sensi- 
tive to  the  assumed  dependence  of  electron  collision  frequency  on  the 


electron  velocity.  The  Appleton-Har tree  equations  assume  the  collision 
frequency  to  be  independent  of  electron  energy  while  the  Sen-Wyller 
equations  assume  the  collision  frequency  to  be  directly  proportional  to 
electron  energy.  Adams  [l97l]  has  shown  that  neither  theory  is  accurate 
below  an  altitude  where  the  collision  frequency  is  greater  than  1.2 
times  the  wave  frequency.  This  limitation  was  based  on  the  comparison  of 
electron  density  profiles  deduced  from  Faraday  rotation  measurements  and 
from  model  calculations. 

The  Faraday  rotation  for  the  wave  propagation  illustrated  in 
Figures  9 and  12  is  plotted  in  Figure  13.  The  altitude  where  = 1.2f 
is  denoted  and  occurs  at  67.2  km.  The  altitude  that  corresponds  to 
a 1“  resolution  in  the  Faraday  rotation  is  well  below  67.2  km.  This 
means  that  Faraday  rotation  data  could  be  distinguished  to  lower  alti- 
tudes, but  unfortunately  interpretation  of  tlie  data  is  not  straightfor- 
ward. For  the  7 MHz  signal  under  consideration,  the  rotation  data  is 
only  good  between  67.2  km  and  72.9  km. 

A discussion  of  the  returns  in  data  by  increasing  the  transmitted 

power  is  now  in  order.  Consider  the  preceding  example,  wliere  it  has 

been  determined  that  differential  absorption  data  down  to  about  46  km 

can  be  secured.  Multiplying  the  transmitted  power  1’.^,  by  a factor  of  10 

to  2 Kw  would  have  the  following  effects.  From  (35)  an  extra  15  dB  range 

would  be  necessary  since  has  been  multiplied  by  a factor  of  about 

4 (see  (34)  ).  The  minimum  signal  drops  off  quickly  at  these  low 

altitudes.  Tlie  effect  on  h^^  determined  by  (33)  would  be  equivalent  to 

lowering  T,  bv  an  order  of  magnitude  (see  Figure  12).  Thus  h would  be 
A ' m 

al>out  2 km  lower.  Barring  the  influence  of  large  man-made  radio  fre- 
quency transmissions,  there  wouldn't  be  any  gains  in  differential  absorp- 
tion information  at  higher  altitudes,  or  Faraday  rotation  throughout  the 
med ium. 

At  the  altitudes  where  both  Faraday  rotation  and  differential  absorp- 
tion data  are  available,  it  is  possible  to  determine  electron  density  and 
the  electron-neutral  collision  rate.  Below  this  common  altitude  region, 
differential  absorption  data  are  still  available  but  independent  electron 
density  and  collision  frequency  determinations  are  not  possible.  The 
collision  frequency  is  generally  thought  to  be  the  most  predictable 


41 


quantity  and  so  is  often  extrapolated  to  lower  altitudes  in  order  to 
determine  electron  densities.  The  procedure  of  extrapolating  tiie  col- 
lision rate  to  lower  altitudes  by  assuming  a particular  pressure  varia- 
tion with  altitude  Pack  and  Phelps,  [l96i]  is  described  in  detail  by 
Thnine  and  Piggott,  [1966]  and  Diakenson  ct  al.,  [l976]. 

The  remainder  of  this  chapter  will  deal  with  the  selection  of  a 
set  of  frequencies  to  probe  the  D-region  under  the  extreme  conditions 
that  are  found  in  the  polar  zone.  Sets  of  curves  corresponding  to  the 
values  of  electron  density  and  collision  frequency  shown  in  Figures  7 
and  8 have  been  developed.  A set  of  curves  of  differential  absorption 
and  Faraday  rotation  for  frequencies  between  3.25  and  25  MHz  have  been 
calculated  based  on  the  disturbed  D-region  electron  density  and  winter 
profiles.  These  are  plotted  in  Figures  14  and  15.  Frequencies  below 
3.25  MHz  were  not  considered  in  this  report  because  it  was  felt  that 
antenna  construction  would  become  too  costly.  One  of  the  goals  of  this 
study  was  to  develop  an  experiment  which  would  be  cheap  and  portable, 
and  hence  would  not  require  great  towers  for  the  transmitter  antenna 
(see  Chapter  IV),  The  3.25  MHz  frequency  should  not  he  considered  a 
limitation  since  our  computer  code  can  generate  Faraday  rotation  and 
differential  absorption  curves  for  any  frequency  of  interest.  Noise 
depends  on  the  time  of  day,  weather,  season  of  year,  and  geographical 
location.  As  a compromise,  the  curve  called  "test  noise"  in  Figure  11 
was  used  to  compute  the  lower  altitude  limit  for  differential  absorption 
measurement . 

The  magnitude  of  differential  absorption  steadily  increases  with 
increasing  altitude  in  Figure  14.  This  corresponds  to  the  absorption  of 
the  extraordinary  mode.  The  absorption  quantities  actually  determined 
from  the  raw  radio  wave  data  are  the  maximum  and  minimum  signal  strengths, 
or  the  modulation  envelope.  The  rate  of  change  in  the  signal  strength 
data,  from  which  dA/dz  is  indirectly  obtained,  is  dependent  upon  the  rate 
of  attervatlon  of  the  extraordinary  mode.  This  can  be  visualized  as  the 
rate  of  change  of  the  signal  envelope  relative  to  the  maximum  signal  (see 
Figures  9 and  10).  Although  dF/dz  hecomevs  larger  with  altitude,  the 
decreasing  change  in  the  signal  envelope  width  relative  to  the  maximum 
signal  implies  that  the  values  of  dA/dz  decrease  in  accuracy  witli 


altitude.  This  problem  is  further  complicated  by  additional  modulation 
of  the  signal  envelope  produced  by  rocket  coning  as  described  in  Appen- 
ilix  C. 

The  upper  altitude  limit  for  differential  absorption  was  found  to 
be  -20.2  dB  and  intersects  the  9 Mllz  curve  at  about  86  km.  A frequency 
slightly  greater  than  9 MHz  remains  within  this  limit  because  of  the 
asymptotic  character  of  the  curves.  As  anticipated,  the  3.25  MHz  sig- 
nal has  the  lowest  altitude  of  measurability  (=46  km)  but  only  propagates 
to  about  66  km  before  the  extraordinary  mode  absorption  is  so  great  that 
the  signal  modulation  can  no  longer  be  measured. 

As  the  rocket  moves  into  the  upper  D-region,  the  rate  at  which  the 
recorded  signal  varies  with  altitude  will  give  an  indication  of  the  reso- 
lution in  the  reduced  data.  It  is  desirable  to  specify  analytically  the 
decrease  in  accuracy  of  differential  absorption  measurements  with  alti- 
tude. We  define  Q as  the  rate  of  change  in  the  signal  envelope  width 
relative  to  the  maximum  signal 


Q = 2 + E ) 

dz  / O X 


I’his  quantity  is  directly  related  to  the  accuracy  of  the  measurement 
and  equivalently,  tlie  number  of  data  points  that  can  be  determined  over 
an  altitude  z.  Ihus,  if  <1  for  a given  freiiuency  is  about  0.1  over  some 
distance,  there  would  be  roughly  ten  times  the  number  of  data  points  as 
comp.ired  to  a signal  with  Q = 0.01  over  the  same  distance.  The  dashed  Q 
contour  lines  plotted  in  Figure  14  indicate  the  relative  accuracy  of  dif- 
ferenti.il  absorption  data  for  various  frequencies.  For  the  D-region 
conditions  represented  by  this  figure,  a transmitter  frequency  of  around 
LO  or  11  MHz  will  give  the  highest  values  of  Q to  the  highest  altitude. 

It  is  seen  from  the  dashed  contours  that  Q decreases  with  increasing 
.iltitude,  and  that  for  any  given  altitude  there  exists  a distinct  fre- 
quency giving  the  largest  value  of  Q.  This  will  be  the  frequency  that 
provides  the  best  resolution  and  hence  the  greatest  accuracy. 

The  corresponding  Faraday  rotation  values  arc  shown  in  Figure  15. 

The  lowest  height  for  measurement  of  Faraday  rotation,  2^,  decreases 
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as  frequency  increases.  The  1°  resolution  forms  a lower  limit,  however, 
for  this  decrease.  At  18-19  MHz  these  two  limits  occur  at  about  59  km. 

The  rate  of  change  of  Faraday  rotation  with  altitude  dF/dz  is  the  quantity 
extracted  from  the  telemetry  data.  Although  dF/dz  is  not  particularly 
small,  it  is  seen  that  the  values  do  indeed  decrease  with  increasing  fre- 
quency. Although  the  differential  absorption  becomes  very  difficult  to 
read  in  the  upper  D-region,  the  electron  density  is  sufficiently  high  that 
the  Faraday  rotation  is  easily  measured.  This  is  the  reason  the  Faraday 
rotation  can  be  recorded  more  accurately  at  higher  altitudes  than  the 
differential  absorption.  Examination  of  the  3.25  MHz  curves  in  Figures  14 
and  15  shows  that  the  upper  limit  for  differential  absorption  measurements 
is  very  near  the  lower  altitude  limit  for  Faraday  rotation.  Thus,  at 
3.25  MHz  there  is  not  altitude  where  simultaneous  measurements  of  Faraday 
rotation  and  differential  absorption  can  be  made. 

Enough  limits  have  now  been  prescribed  that  a choice  of  frequencies 
can  be  made  to  achieve  optimum  results  for  expected  ionospheric  condi- 
tions. For  a disturbed  polar  ionospliere  and  moderate  collision  rates 
such  as  found  in  winter,  the  following  frequencies  are  recommended. 

19.0  MHz  - Differential  absorption  data  for 

48  to  85  km.  Faraday  rotation  data 
for  59  to  110  km. 

10  MHz  - Differential  absorption  data  for  48 
to  90  km.  Faraday  rotation  for  65 
to  110  km. 

3.25  MHz  - Differential  absorption  data  for  46 
to  66  km.  No  Faraday  rotation  data. 

The  upper  altitude  limit  for  differential  absorption  was  chos-. 
for  a value  of  Q which  would  give  one  data  point  per  kilomi'ier.  It 
should  be  emphasized  that  above  this  level  there  is  still  d if f t'rent iai 
absorption  that  can  be  measured,  but  the  rate  of  absorpti(’n  is  too  low 
to  yield  an  accurate  value  every  kilometer. 

It  is  not  possible  to  determine  both  electron  density  and  collision 
frequency  below  59  km  since  that  is  the  lowest  altitude  for  Faraday  rota- 
tion measurements.  The  collision  frequency  can  be  extrapolated  downward 
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to  obtain  estimated  electrcm  densities  below  59  km.  Even  though  it  is 
indicated  that  Faraday  rotation  data  for  19  MHz  can  be  read  down  to 
48  kin,  the  data  is  a factor  of  2 to  3 less  accurate  than  the  10  MHz 
signal  below  70  km. 

Receiver  requirements  for  this  experiment  are  greatest  for  the 
3.25  MHz  signal  where  the  noise  temperature  T = 3 x 10^°K  corresponds 
to  0.38  pV/m.  The  lowest  range  full  scale  voltage  E must  then  be 
about  4.2  pV/m.  At  h^  = 46  km,  E^  =6.8  mV/m,  so  that  4 receiver  ranges 
are  needed.  (As  noise  increases,  less  sensitivity  is  useable  in  the 
receiver . ) 

The  calculations  described  above  were  repeated  for  a disturbed 
D-region  but  with  the  high  collision  frequency  profile  of  Figure  8. 

The  values  are  plotted  in  Figures  16  and  17.  The  lilgher  collision  rate 
increases  the  total  absorption  with  the  net  result  that  the  lowest  alti- 
tude for  reliable  measurements  has  been  raised  about  4 km.  For  the  same 
reason  higher  frequencies  must  be  used  to  measure  the  top  portion  of  the 
l)-region  ionosphere.  The  dashed  lines  in  Figure  16  represent  constant 
values  of  Q.  The  1°  and  limits  for  Faraday  rotation  measurements 
occur  about  65  km  for  a 25  MHz  wave,  Figure  17. 

Tile  three  optimum  frequencies  for  this  model  are: 

25  MHz  - Differential  absorption  data  for  54  to 
90  km.  Faraday  rotation  data  for  65  to 
no  km. 

14  MHz  - Differential  absorption  data  for  51  to 
96  km.  Rotation  data  for  71  to  110  km. 

3.25  MHz  - Differential  absorption  data  from  about 
48  to  70  km. 

Receiver  requirements  for  these  conditions  are  similar  to  the 
disturbed  winter  case.  For  an  equivalent  electric  field  noise  of 
0.38  pV/m,  4 receiver  ranges  are  needed  to  measure  the  signal  at  51  km 
where  the  minimum  signal  is  1.2  pV/m  and  the  maximum  signal  is  6.2  mV/m. 

The  third  case  oi  interest  is  that  of  a quiet  D-region  and  winter 
collision  rates.  Utilizing  the appropr  ia  t e curves  in  Figures  7 and  8, 
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DIFFERENTIAL  ABSORPTION  (dB) 


t'igiiro  1().  I)i  fforcMU  ial  absorption  profiles  for  the  disturbed 
electron  densitv  and  summer  collision  frequencv 
profiles  of  Figures  7 and  8. 
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FARADAY  ROTATION  (Degrees) 


ri)',iiri'  17.  Knr.idiiy  rot.jtJon  profilis  for  the  disturbed 

olertron  density  and  svimner  collision  frequency 
protiles  of  Figures  7 and  8. 
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the  differential  absorption  and  Faraday  rotation  were  calculated  and 
plotted  in  Figures  18  and  19.  Values  of  dA/dz  and  dF/dz  for  this  case 
are  much  larger  than  for  the  disturbed  cases.  For  the  case  of  Faraday 
rotation,  z^  for  all  frequencies  is  well  below  the  l”  Faraday  rotation 
altitude.  For  differential  absorption  measurements  to  110  km,  a fre- 
quency close  to  4 MHz  is  needed.  Frequencies  greater  than  7 MHz  have 
values  Ilf  Q less  than  0.05  over  the  entire  range.  For  this  reason,  they 
are  not  included  in  the  figures.  For  measurements  under  these  ionosi)lier  ic 
conditions  two  frequencies  would  be  selected; 

4.25  MHz  - Differential  absorption  data  for 

75  to  110  kill.  Faraday  rotation 
data  81  to  110  km. 

3.25  MHz  - Differential  absorption  data  for 

73  to  103  km.  Faraday  rotation 
data  for  79  to  106  km. 

A 2.5  MHz  curve  was  calculated  for  comparison  with  experiments  described 
by  Pi'HHi'tt  L'i  [l972]  and  Li: , [1974]  under  similar  atmospheric 

conditions.  As  can  be  seen  this  lower  frequency  yields  more  accvirate 
results  over  the  70-90  km  height  region.  Ri-ceiver  requirements  for  t lie 

3.25  MHz  signal  at  the  lower  limit  of  73  km  also  call  for  4 ranges,  tlie 
highest  me.tsuring  a field  of  4.4  millivolts/meter.  This  gives  a total 
dynamic  amplitude  range  of  2811:1. 

In  conclusion,  it  can  be  seen  by  comparison  of  Figures  14  and  18 
and  Figure  7,  that  while  electron  density  measurements  as  low  as  10  cm~^ 
are  possible  under  quiet  conditions,  barely  100  cm“^  are  attainable  for 
the  disturbed  case.  The  difference,  of  course,  is  the  higher  degree  of 
absorption  for  the  disturbed  case  for  any  given  electron  density.  This 
dependence  also  shows  up  by  comparing  Figures  14  and  16.  Thus,  as  a 
general  rule,  measurement  accuracy  and  electron  density  determination 
at  lower  altitudes  decrease  with  increasing  electron  collision  frequency 
prof i le. 


CHAPTER  IV 


S3 


ANTENNAS 

Ground  Station  Antennas 

Signals  from  the  ground  station  must  radiate  throughout  that  por- 
tion of  the  ionosphere  traversed  by  the  missile.  A radiation  pattern 
similar  to  that  of  a one-half  wavelength  dipole  mounted  horizontally 
one-fourth  wavelength  above  a ground  plane  would  be  suitable.  The 
antennas  must  be  easily  erected  and  movable  from  location  to  location 
as  needed.  It  is  to  be  noted  that  supporting  masts,  to  be  practical, 
cannot  exceed  about  25  or  30  feet  in  lieight  without  being  permanently 
placed  in  position.  The  minimum  frequency  that  could  easily  be  used 
(with  A/4  height)  would  then  be  about  8 MHz.  However,  horizontal  di- 
poles may  be  mounted  somewhat  lower  than  one- fourth  wavelengtii  without 
seriously  affecting  the  radiation  pattern. 

Kraus  [l950]  gives  a relationship  for  the  gain  in  field  intensity 
at  large  distances  for  a "Half-Wavelength  Antenna  Above  Ground"  (HWAG) 
with  respect  to  that  of  a "Half-Wavelength  Antenna  in  Free  Space"  (HWFS) 


as 
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/' 
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_ 
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h^  = 2tt/A  times  the  height  above  ground  = (2'7X)h 
Rj^j  = self  resistance  of  a >/2  antenna 


11. 


loss  resistance  of  the  'V2  antennas 


R = mutual  resistance  of  the  A/2  antenna  and  its  image  at  a 
distance  of  2 h 

Plots  of  radiation  patterns  are  given  in  Figure  20  for  li  = 0.25  X,  0.1  A, 
and  0.05  X.  R^^^  is  assumed  zeros  for  the  plots  shown.  It  is  seen  that 
as  h becomes  smaller,  the  field  intensity  pattern  deviates  only  slightlv 
from  the  pattern  for  h = A/4.  However,  the  terminal  impedance  decreases 
rapidly  as  h decreases.  For  plot  A the  terminal  impedance  is 
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Horizontal  half 
wave  dipole 


Ground  (Ideal) 


The  half-wavelength  dipoles, 
plane  as  in  (a)  produce  field 
shown  in  (b)  for  A,  h = A/4; 


mounted  above  a ground 
intensitv  patterns  as 
B,  h = A/ 10:  and  (1, 


approximately  86  ohms,  for  B it  is  about  22  ohms,  and  for  C it  is  about 
6 olims.  Wiien  h is  very  small,  the  effieiency  of  the  antenna  would  de- 
crease because  the  loss  resistance  becomes  significant  compared  to  the 
radiation  resistance.  Witli  minimum  h chosen  as  X/10,  and  with  masts 
no  iiigher  than  30  feet,  it  is  seen  that  the  minimum  frequency  to  be 
used  for  the  propagation  experiment  should  not  be  lower  than  about  3.25  MHz. 
As  the  frequency  is  increased,  the  decreasing  wavelength  will  permit  the 
antenna  to  be  mounted  a greater  fraction  of  a wavelength  above  the  ground 
p 1 ane . 

Matcliing  the  antennas  to  the  coaxial  cable  feed  lines  may  be  accom- 
plished by  using  either  of  several  imitching  networks.  The  lowest  fre- 
quency antenna  (h  = A/10)  with  a terminal  impedance  of  = 22  + jx  for 
a A/2  antenna  may  be  efficiently  matched  rather  easily.  mav  be 

increased  by  a factor  of  9 by  using  a 3-wire  folded  dijiole.  As  shown  in 
Figure  21,  a coa.\ial  cable  drives  one-half  of  the  antenna  directly  from 
the  center  conductor.  The  otlier  half  of  tiie  antenna  is  excited  from  the 
center  conductor  through  a A/2  length  of  coaxial  cable.  The  resulting 
•erminal  impedance  is  approximately  49  + jO  ohms  after  the  length  of  the 
antenna  lias  been  adjusted  for  zero  reactance. 

Wlien  higher  frequencies  are  cliosen,  the  height  of  the  antennas  may 
be  adjusted  to  li  ’ A/10,  and  then  2-wire  folded-dipoles  may  be  used.  The 
terminal  impedance  may  again  be  adjusted  for  an  excellent  match  with  the 
50  ohm  source  cable,  as  in  Figure?!,  by  adjusting  the  height  above  the 
ground  plane  and  by  trimming  the  antenna  length  for  zero  reactance.  The 
height  of  the  antenna  will  then  be  about  0.16  A above  a well  defined 
ground  plane. 

Alternate  matching  networks  make  use  of  the  balun  for  matching 
balanced  to  unbalanced  impedances,  or  the  network  in  Figure?’ may  be  used. 

Kach  of  the  matching  methods  described  for  the  folded-  and  single- 
wire antennas  have  been  built  and  tested  successfully  at  a frequency  of 
4 MHz  in  connection  with  this  report. 

The  theoretical  field  intensity  patterns  shown  in  Figure  20  are 
difficult  to  measure  without  a suitable  antenna  range.  However,  an 
attempt  was  made  to  verify  the  radiation  patterns  for  h = A/4  and  A/10, 
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3-wlre  folded  X/2  antenna 


Figure  21.  A three-wire  folded  half-wavelength  antenna  phased 
and  driven  by  a coaxial  cable.  With  h = A/10  above 
the  ground  plane  and  with  the  length  trimmed  for 
zero  reactance,  the  terminal  impedance  presented  to 
the  transmission  line  is  Z^.  = 22  x (3)^/4  = 49  Q. 


Figure  22.  A tuned  matching  network  for  matching  an  unbalanced 
signal  source  to  a balanced  dipole  antenna. 
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Vertical  dipoles  were  mounted  in  front  of  a vertical  screen  and  the  radi- 
ation in  the  ground  plane  was  measured.  Field  intensity  patterns  wore 
measured  with  the  spacing  between  the  antenna  and  the  screen  adjusted 
first  to  X/10  and  then  to  A/4.  The  model  was  operated  at  24  MHz  which 
permitted  use  of  tiie  smallest  possible  model  within  the  range  of  avail- 
able field-intensity  measuring  equipment.  The  vertical  screen  was  40  feet 
long  and  15  feet  high  and  consisted  of  vertical  wires  spaced  1 foot  apart. 
Since  the  screen  was  only  about  1 wavelength  long,  the  horizontal  radia- 
tion pattern  was  not  expected  to  match  the  ideal  pattern  but  only  to 
verify  the  performance.  The  antenna  and  screen  configurations  are  shown 
in  Figure  23.  The  measured  normalized  field  intensity  patterns  are  shown 
in  Figure  24.  The  field  intensity  was  measured  45  meters  from  the  center 
of  the  screen. 


Figure  23.  Tiie  field  intensity  measurements  were  made  in  the 

ground  plane  around  the  dipole  and  the  screen.  Tiie 
antenna  was  excited  through  a coaxial  cable  from  a 
signal  generator  located  at  x behind  the  screen. 
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An  t ejinas_  JAi  r_  t lie  borne  Keco  i ve  r 

The  small  size  and  high  velocity  of  th€=  sounding  rocket  make  the  use 
of  outboard  antennas  impractical.  llierefore,  ferrite  core  antennas  are 
proposed  for  use  with  the  payload. 

Ferrite  core  antennas  were  fabricated  and  tested  at  frequencies  of 
4 MHz,  17  MHz,  and  25  MHz.  Initial  tests  wore  made  using  ferrite  rods 
measuring  0.337  inches  diarater  x 7.5  inches  long  of  and  material 
from  Permag  Sierra  Corporation.  Tlie  objective  was  to  determine  the  vol- 
tage at  tile  output  terminals  of  the  antennas  versus  the  field  intensity. 
'Ilie  Qj  material  was  clearly  superior  in  performance  at  4 MHz,  the 
material  gave  tlie  best  results  at  25  MHz,  while  at  17  Mliz  the  two  ma- 
terials gave  approximately  equal  results. 

The  final  tests  were  made  witli  the  same  core  material  shortened  to 
3.75  inches.  The  results  using  the  shortened  ferrite  rods  were  approxi- 
mately the  same  as  those  observed  witli  the  longer  antenna  rods  except 
that  the  terminal  voltages  were  rouglily  one-third  less.  It  is  antici- 
pated that  antenna  rods  will  be  less  than  4 inches  in  length  for  most 
payloads  using  the  propagation  experiment  because  of  the  likelihood  that 
small  missiles  will  be  used. 

The  antennas  were  fabricated  by  winding  the  ferrite  rods  with  closelv 
spaced  turns  of  No.  24  Nylclad  copper  wire.  Each  antenna  was  tuned  with 
about  iO  pK  capacitance.  Tlie  terminal  impedance  of  the  antennas  was  ad- 
justed to  1,000  ohms.  An  emitter  follower  was  then  used  to  match  the 
antenna  impedance  to  the  50  ohm  terminal  impedance  of  a Stoddard  Model 
No.  NM20-B  RI-FI  Meter  whicli  served  as  the  field  intensity  meter  and 
also  as  a sensitive  voltmeter.  Figure  25  shows  the  experimental  arrange- 
ment used  for  measuring  antenna  terminal  voltage  versus  field  intensity. 

A typical  result  at  4 MHz  was  65  iiV  delivered  by  the  antenna  to  a 
1,000  ohm  load  with  1.55  mV/meter  of  electric  field  intensity.  It  is 
safe  to  say  that  1 pV  of  signal  would  produce  a favorable  signal  to 
noise  ratio  so  that  a signal  strength  of  1.55  x 10“^  (V/m)/65  = 24  x 
10“'’  V/m.  If  the  greatest  distance  between  the  transmitting  and  receiving 
antennas  is  10’  meters  and  if  the  gain  of  the  transmitting  antenna  is 
given  a safe  value  of  3.5,  then  we  may  calculate  the  iransmitted  power 


60 


I'igiiro  25.  The  experimental  test  used  to  determine  antenna 
terminal  voltage  versus  field  intensity.  Both 
the  output  voltage  of  the  antenna  and  the  strength 
of  the  E-field  were  measured  witli  the  IR-Fl  Meter. 


without  the  effects  of  Faraday  rotation  and  attenuation  through  the  iono- 
sphere. There  is  good  evidence  from  the  work  of  Dr.  Harris  (Utah  State 
University)  that  tlie  ilynamic  range  imposed  by  Faraday  rotation  is  about 
100  to  1,  and  that  the  absorption  of  the  signal  at  the  lowest  frequencies 
of  interest  imposes  another  factor  of  100.  Tlie  transmittal  power  may  then 
be  calculated  from  the  relat  ionsliip 

4Tir ^ P j 

~ ^ (attenuation  factor)  x (Faraday  rotation  factor) 


P - power  density  at  the  receiving  antenna 


E‘  ^ (24  X 10  V_/m)^ 

2.  'W?  ' 

o 


1.5  3 X 10  W/m^ 


1 


r = 1 0 ^ m 


niat  level  of  transmitted  power  Is  not  only  reasonable  but  it  appears  to 
be  on  the  safe  side  since  a good  receiver  can  operate  with  a suitable 
signal-to-noise  ratio  with  about  0.1  mV  of  input  signal  in  a relatively 
quiet  environment.  Since  the  environmental  noise  at  the  receiver  is  a 
composite  from  many  possible  sources  and  is  unknown  a priori,  it  would 
seem  that  the  ground  based  transmittei  should  deliver  700  to  1,000  watts 
to  the  transmitting  antennas. 

A^l  though  the  transmitted  power  requirement  may  vary  slight  Iv  from 
one  frequency  to  another,  the  worst  case  is  that  of  the  example  above. 
Therefore,  transmitted  power  in  the  range  of  700  to  1,000  watts  should 
be  adequate  for  any  frequency  utilized  in  the  propagation  experiment. 
Inexpensive  transmitters  are  readily  available  for  these  power  levels 
and  the  frequencies  of  interest. 
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CHAPTER  V 
CONCLUSIONS 

The  goal  of  this  work  was  the  development  of  a multiple  frequency 
propagation  experiment  to  measure  electron  densities  and  electron  colli- 
sion rates  in  the  ionospheric  D-region.  This  experiment  was  designed  to 
utilize  standard  HF  transmitters  and  simple  antenna  systems  that  could 
be  low  power,  portable  and  inexpensive.  This  experiment  employs  the 
radio  wave  signals  to  measure  the  Faraday  rotation  and  differential 
absorption  between  ground  transmitters  and  rocket-borne  receivers  tra- 
veling through  the  D-region.  The  time  history  of  the  Faraday  rotation 
and  differential  absorption  then  permit  an  evaluation  of  the  height 
profiles  of  electron  density  and  collision  rate. 

A theoretical  propagation  model,  based  on  the  Sen-Wyller  form  of 
the  magnetoionic  equations,  was  developed  to  numerically  calculate  the 
wave  propagation  from  ground  to  rocket  at  high  magnetic  latitudes. 

This  propagation  model  included  the  limitations  imposed  by  noise,  equip- 
ment sensitivity  and  data  evaluation.  Curves  of  Faraday  rotation  and 
differential  absorption,  based  on  extreme  values  of  electron  density 
and  electron  collision  rate  as  might  be  seen  in  the  polar  lat  itudes , were 
calculated  to  illustrate  the  frequencies  and  altitudes  at  which  this 
propagation  experiment  would  be  successful.  It  was  found  that  under 
quiet  ionospheric  conditions  electron  densities  as  low  as  10  cm“^  could 
be  measured  by  this  propagation  experiment,  whereas  increased  collision 
rates  such  as  occur  in  disturbed  conditions  reduce  tlie  resolution  to 
less  than  100  cm~^.  Appendix  C contains  an  extensive  derivation  of  the 
effect  of  rocket  coning  on  the  received  signal.  We  found  that  rocket 
coning  could  produce  significant  errors  in  the  apparent  values  of  Faraday 
rotation  and  differential  absorption,  and  that  these  errors  were  not 
systematic  and  not  removable  by  simple  averaging  technique.  A metitod 
was  developed  to  remove  this  perturbation  from  the  Faraday  rotation  and 
differential  absorption  data  with  the  aid  of  continuous  rocket  attitude 
information.  This  technique  also  permits  the  calculation  of  electron 
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density  and  electron  collision  rate  from  the  raw  Faraday  rotation  and 
differential  absorption  data. 

Chapter  IV  summarizes  a study  of  antennas  proposed  for  use  in  the 
rocket  and  also  for  the  ground  transmitter.  It  appears  that  a trans- 
mitter power  in  the  range  of  700  to  1000  watts  radiating  on  horizontal 
half-wave  dipoles  will  be  adequate  for  this  propagation  experiment.  The 
rocketborne  receiving  antennas  utilize  ferrite  core  loops. 

The  calculated  curves  of  Faraday  rotation  and  differential  absorp- 
tion together  with  the  constraints  already  mentioned  were  used  to  evalu- 
ate wave  frequencies  for  the  extreme  ionospheric  conditions  posed  in 
Chapter  III.  These  are  summarized  in  the  following  table.  Although 
1.25  MHz  was  the  lowest  frequency  considered  in  this  study,  frequencies 
lower  than  3.25  MHz  miglit  be  useful  in  other  special  cases. 


I’olar  D-Kegion 
Ionospheric  Conditions 

Proposed 

Frequencies  (MHz) 

Winter  disturbed 

3.25,  10.,  19. 

Winter  quiet 

3.25,  4.25 

Summer  disturbed 

3.25,  14.,  25. 

Summer  quiet 

3.25,  6.0 

It  is  our  conclusion  that  a mul t i f requency  prc'pagation  experiment  is 
a viable  scheme  for  obtaining  accurate  D-region  electron  density  profiles. 
The  techniques  that  have  been  developed  in  this  work  allow  an  evaluation 
ol  frequencies  before  hand  so  tliat  optimum  results  can  be  obtained;  and 
accurate  nu-thods  of  data  evaluation  after  the  rocket  flight. 
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APPENDIX  A.  MAGNETOIONIC  FORMUEAS 


Tho  generalized  complex  refractive  index  given  hy  and  Wulh-r', 
[i960]  is 


('ll  _ _ A + B stn^  +v/B^sin‘''l)  - c^cos^ii 

D + E sin^'d) 


where 


A=2.i(>^.il^l) 

" ’ I 1 1 I *^111^  ''‘'ii 

"=2^1-11 
I)  = Tr 


E = 2r 


Ihe  f's  are  elements  of  the  generalized  dielectric 
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, j = '^[(f-d)  -t  l(c-e)  ] 

■ .,=  - '^(c+e)  + i[b  - !i(f+d)] 


a = -^7  C , 
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■V? 


= V 0 

2..,v„  >/-- 


Li)  + S 
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liu-  clualLtios  f>',  UV,,  are  wave  freqiioncv,  plasma  frequency,  and 
collision  frequency,  respectively.  The  semiconductor  integrals 
represented  by  C3/3  and  are  described  by  Di’mlp,  of  ,7/.  [l9S7] 

and  can  be  suitably  approximated  by  polynomials  given  by  Unva  [196T], 
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4 

X 

+ b,v 


+ <1? 
+ b,,  x'*  + b 


+ a 
’ + b 




+ h,x  + b„ 


•Tj,  = 2.3981474  x 10'^ 

Uj  = 1.1287513  X 10 

a^  = 1,1  394160  x 10"^^ 

a,  = 2.46531 15x10 


h,  = 1 .8064128  X 10"^ 
b,  = 9.3877372 
bj  = 1.4921254  X 10+2 
hj  = 2.8958085  x 10+^ 
= 1.204951  2 X lO"^'" 
85  = 1,4656819  X 10 


C,//x) 


+ a + a, 

b,.x'*  + b,x-^  + bjx^  + b,x  + bp 


1 . 16  5064! 
1.6901002  X 10 
6.6945939 


bj,  = 4.3605732 
6j  = 6.4093464  x 10 
= 6.8920505  x 10 
63  = 3.5355257  X 10 
b^  = 6.6314497 


APPENDIX  B 

THE  WAVE  ELECTRIC  El  ELD  VECTOR 


Accurate  reduction  of  telemetered  propagation  data  requires  a deter- 
mination of  the  ray  path  between  the  radio  transmitter  and  the  rocket 
receiver.  in  gem-ral,  an  electromagnetic  wave  propagating  through  a con- 
tinuously varying,  stratified  anisotropic  medium  will  bend  away  from 
regions  where-  tlie  index  of  refraction  \i  decreases.  This  phenomenon  leads 
to  the  process  oi  radio  wave  reflection  from  the  ionosphere.  (See 
pp.  190-191). 

In  order  to  caliulate  the  ray  path  it  is  general Iv  necessary  to 
liivide  the  ion  i tted  medium  into  thin  layers  with  the  requirement  that 
~ D-,  “1*.  and  .Ay.  = X-.  "X-  't''*-’  small  compared  to  y..  The 
i ^ layer  may  then  be  ch.iract  eritted  by  average  vtilues  of  the  indices 
. = (D-.  + 9.)/h  iind  X.*=  (X.j  + ^X-,  + y.)/fi,  where 

, and  / .ir-'  the  indices  at  (z,  + z . , )/2.  The  change  in  rav 

patli  direction  at  the  surface,  of  two  media  of  different  refractive  in- 
dex is  given  bv  Snell's  law;  y.sinO.  = ii . sin  , where  the  variation 

1 1 i+i  1+1 

in  V.'  over  the  r.iy  p.ith  is  small  enough  (less  than  5°)  to  be  ignored  so 

that  D.,  ( - e.,  (i^.)  • 

1+1  1+1  1+1  1 

Utilizing  this  expression  to  calculate  the  deflection  of  a typical 
w.ivi'  frequency  propagating  through  the  altitudes  where  d i f f er  .-n  t ia  1 
absorption  c.in  be  mt-asured,  it  w.is  found  that  in  a A km  thick  layer, 
tile  deilei'tion  was  less  than  0.1  . Thus,  in  our  c.il  cu  la  t ions  of  dif- 
ferential .absorption  and  Earaday  nit.ation  we  can  assume  tliat  the  rav 
p.itlis  are  str.iight  lines. 

Anotlier  char.ictor  ist  ic  of  propagating  waves  in  an  ionized  medium 
is  the  tendency  for  the  E-field  to  acquire  a component  in  the  direction 
of  propagation.  For  c.alculations  involving  rocket  coning,  it  is  neces- 
sary to  show  that  the-  longitudinal  Ji-field  component  is  very  small.  To 
determine  the  plane  of  the  phasor  _E,  we  must  find  the  relative  magnitude 
t>f  E^  (the  component  of  E in  the  direction  of  wave  propagation)  in  terms 
<’f  E^  or  E^.  This  can  be  found  by  a derivation  similar  to  that  for  the 
wave  polarization  R.  By  substituting  the  dielectric  tensor  in  Appendix  A 
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into  M.ixwt'll's  equations  and  solving  for  E /E  , we  can  numerically 

z y 

determine  the  relative  magnitude  of  the  longitudinal  E-field  component 
for  the  coordinate  system  of  Figure  3.  The  result 


'z  ^ [(M  - jX)^  E - 2b]  sint}) 


can  be  evaluated  hv  substituting  tvpical  values  for  w,  v,  n , etc.  into 

e 

the  expressions  for  B,  C,  E,  and  (^.  Calculations  show  that  E /R  < 

z y “ 

10  + ) 10  for  typical  ionospheric  parameters. 

We  can  safely  conclude  that  for  all  practical  applications  E lies 
in  a plane  perpendicular  to  the  ray  path,  which  was  shown  to  be  essen- 
tially a straight  line  from  the  ground  antenna  to  the  rocket  vehicle. 
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The  effects  of  coning  on  the  signal  at  the  rocket  can  be  described 
analytically  by  developing  equations  tliat  transform  the  coordinate  frame 
in  which  the  antenna  is  stationary  to  that  in  which  the  ellipse  traced 
out  by  the  radio  wave  E-field  vector  is  stationary,  the  Faraday  rotation 
frame.  Coordinate  rotations  are  required  that  transform  (x^ , , z'^) , 

where  z‘*  is  defined  as  the  rocket  spin  axis,  is  the  direction  along 
the  antenna,  and  y‘^  is  perpendicular  to  both  ^ and  to  the  coordinates 

(x^- , , where  is  the  coning  axis,  y*-  points  magnetic  north  and 

y*-  is  perpendicular  to  bothy*-  and  . The  rotational  motion  of  the 
rocket  in  a frame  fixed  relative  to  an  observer  on  the  earth,  can  be 
describeii  by  a = da/dt,  the  rocket  spin  rate.  The  coning  angle  B is  the 
• ingle  between  the  rocket  major  .ixis  (spin  axis  y'’)  and  the  coning  axis  y*- . 
Kinaily,  define  y = dj/dt  to  be  tlie  coning  rate;  tlie  angular  velocity  by 
which  the  rocket  axis  y'*  processes  aliout  the  axis  z*  . Figure  Cl  shows 
tlie  spatial  relationships  between  the  two  coordinate  frames  and  the  angles 
defined  .ibove.  A clockwise  rotation  ft  about  the  axis  can  be  repre- 
sented by 

( cos  ut  -sin  at  0 ) 

S = sin  at  cos  at  0 (Cl) 

( 0 0 I I 

Thus  = Sy‘* . The  next  rotation  is  about  y*’  and  causes  y*’  = Sy^ 
to  be  rotated  by  an  angle  P into  coincidence  with  y*"  . This  is  a clock- 
wise rotation  that  can  be  represented  by 


' 

(-  = s 0 cos  B 

I 0 sin  B 


-sin  B 
cos  P 


(C2) 


The  geometry  of  this  situation  is  illustrated  in  Figure  C2  where  y*'  = 
Cy*’  = CSy  . The  axis  y'  = Cy^’  = CSy'*  lies  in  the  plane  defined  by  y*" 


and  y 


The  rocket  antenna  now  lies  in  the  (x'  , /)  plane  by  the  trans- 


formations S and  C so  that  x 


Cx^  = CSx‘^ 


Finally,  the  (x^  plane 


Ls  rotatcxi  about  the  z axis  by  y degrees  so  that  x"  and  / now  align  with 
X and  y 

sented  by 


x'  and  respectively.  Again  this  is  a clockwise  rotation,  and  is  repre 


COS  Yt 

P - > • • 

r = Sin  Yt 

0 


-sin  Yt 
cos  Yt 
0 


(C3) 


The  matrix  PCS 
a coordinate'  syste'm 
tile  Vitound. 


transforms  the  stationary  rocket  antenna  axes  into 
that  IS  stationary  with  respect  to  an  observer  on 


Stationary  antenna  frame  (x‘^,  y‘\ 
relative  to  (x'^,  yt.  zti  ~ ’ 


£ ) orientation 


Figure  Cl 


Figure  C2.  Orientation  of  (x  » y , 7.  ) 
to  (x^  yt,  ,.t)  - ^ - 


Assume  .in  .■in>;le  <5  exists  between  llie  r.idio  w.ive  r.iy  paths  and  the 
axis  z .about  whieh  the  rocket  precesses.  The  rotation  matrix 


cos  (S 
-sip  6 


sin 


describes  a rotation  about  the  axis,  taking  the  (y*^,  z*^)  axes  into 
the  (y  , ^ ) .ixes  of  Figure  2 of  Chapter  II.  The  transformation  from 
^21  > y ’ ^ the  larad.ay  rotation  frame  (x'^,  y*^,  was  developed  in 


r 
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Chapter  II  and  is  contained  in  the  matrix  R of  equation  (14).  Summln}’ 
all  the  rotations  gives  the  matrix 

T = RDPCS  (C5) 


which  transforms  the  antenna  coordinate  axes  of  a spinning,  coning  rocket 
to  the  Faraday  rotation  frame  in  the  ionospheric  medium. 

The  angle  6 between  the  ray  path  and  the  coning  axis  may  be 

impossible  to  determine  experimentally.  We  assume  that  points  in  the 
direction  of  the  spin  axis  at  launch.  If  the  transmitter  is  close  to  the 
launch  site,  then  the  assumption  that  6 is  less  than  5°  is  probably  reason- 
able at  least  for  the  up-leg  and  the  D-region.  With  6 less  than  about  30°, 
the  rotation  D can  be  approximated  by  the  expression 


D = 


0 

1 

0 


1 


I 0 0 0 

+ 0 -6‘/2  6 

( 0 -6  -6^/2 


1 + D' 


(C6) 


with  less  than  5Z  error.  Since  the  measurement  of  the  coning  angle 
depends  upon  the  linearity  of  the  magnetometer,  we  can  associate  anv 
non-negl igible  experimental  error  by  b.  Substituting  B + b into  (C2)  with 
b small  allows  us  to  write 


0 0 
ensB  -sinB 

sinB  cosS 


+ 


I 0 0 

I 0 ‘^bsinB 

( 0 ‘bcosB 


‘bsinB  , (C7) 

*bsinB  ) 


~ C,  + E' 


where  E'  is  an  error  matrix.  Substituting  (C7)  and  (C6)  into  (C5)  gives 


T = RPCS  + RD'PCS  + RPE'S  + higher  order  error  terms. 


1 


The  matrix  product  RPCS  is  the  principal  term;  the  other  terms  being  error 
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matrices  wlioso  elements  are  only  a small  percentage  c)f  the  principal 
matrix  elements.  The  matrix  elements  for  T s RPCS  are 

• * 

= cosy  cosat  - siny ' cosPsinat 

♦ • 

T = -cosy  sinat  - siny ' cosgcosat 
1 2 

T = siny'sinB 

1 3 

• • 

T = siny  cosat  + cosy ' cosBsinat 

2 1 

= -siny 'sinat  + cosy ' cosPsinat  (C8) 

T = -cosy'sinB 

2 3 

T = sinPsinat 

3 1 

T = sinScosat 

3 2 

T = cosB 

3 3 


where  y'  = yt  - F and  F = ((*i^  - the  Faraday  rotation  angle. 

(\s  a resvilt  of  ignoring  S,  the  coning  angle  yt  is  measured  in  the 
plane  of  polarization  so  that  the  Faraday  rotation  either  advances  or 
retards  the  precession  effect  appearing  in  the  unreduced  data.  The 
direction  of  the  rocket  dipole  antenna  relative  to  the  Faraday  rotation 
Irame  is  = T • where  is  the  unit  vector  in  the  ^ direction. 

= [cosy 'cosat  - s iny ' cosBsinat 

+ [siny'cosat  + cosy ' cosBs inat]^*^  + sinBsinat  ^ . (C9) 

The  potential  measured  at  the  base  of  the  antenna  is  proportional  to  the 
integrated  electric  field  component  parallel  to  the  antenna.  The  scalar 
product  of  equation  16  of  Chapter  II  and  (C9)  gives 

F. . = [cosy' cosat  - siny ' cospsinatjE  cos[ut  - (<i  + iJi  )/2] 

'ox 

- [siny'cosat  + cosy ' cos8sinat]E2sin[a)t  - ((J)^  + <^^)/2]  (CIO) 


where  F.^  and  E^  are  given  by  equation  19  of  Chapter  II.  In  the  absence 
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of  coning,  = 0 and  y'  combines  with  a 

24  of  Chapter  II.  The  maximum  value  of 

dE  ' 

2TT/a)  is  found  by  setting  — - = 0 with 

dt 

The  result 


so  that  (CIO)  reduces  to  equation 

equation  (CIO)  in  the  wave  period 

• • 

the  assumption  that  ot,  y <<  w. 


wt 


(cj)  + ) 

O X 


E.> 


-tan 


• • I 

/sinY'cosat  + cosy ' eos^sinat  W (cn) 
U'osY'cos(^t  - sinY  ' cos3sinat  / I 


can  be  substituted  back  into  (CIO)  giving 


(<},  + ] 

II 

< 

cos 

0 X 

[cosY^cosc^t  - sinY ' cosBsinctt] 

max 

-f 

< 

s in 

(alt  

[sinY  * cosott  + cosy  * cosSsinat] 

max 

where  = ?-Ej,  V.,  = , and  i is  the  effective  electrical  length. 

Ttto  voltage  in  a coning  rocket  is  equivalent  to  the  voltage  , 
eqiKition  25  of  Chapter  II  in  a non-coning  rocket.  It  is  seen  that  the 
peak  amplitude  of  the  signal  is  modulated  by  the  relative  aspect  of  the 

ray  path  and  spin  axis.  If  the  coning  rate  were  rapid  with  respect  to 

the  rocket  velocity,  then  an  averaging  process  would  smooth  out  coning 
effects.  However,  if  the  altitude  change  in  one  precession  cycle  is  sig- 
nificant, then  Faraday  rotation  and  differential  absorption  also  change 

in  the  i-ycle  and  averaging  will  not  yield  accurate  results.  The  spin 

rate  must  be  determined  experimentally  by  magnetometers  on  the  rocket, 
but  theoretically,  half-rotations  are  solutions  of  the  transcendental 
equa  t ion 


tan  ' t ^ (90°  - yt) 
a L cosP 


(Cl  3) 


and  correspond  to  the  instant  the  antenna  is  parallel  to  the  earth's 
iN.igni’t  ic  tieid.  Equation  (Cl  3)  is  found  by  setting  the  element  of  the 
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first  row  and  lohimn  of  matrix  I’t'Soqu.il  to  zoro.  If  Faradav  rotation 
is  zoro  (>'  = Yf),  t lion  t ho  'spin  marks'  aro  coinoidont  wifli  tho  minimum 
valiios  of  (Cl  2). 

Tho  complexity  of  (Cl  2)  clouds  the  easy  determination  of  Faraday 
rotation  and  ti  i f f eren  t ia  1 absorption.  lYlnit  we  need  then  is  a systematic 
method  of  deducing  tlieso  quantities  in  tlie  ease  of  rocket  coning.  Vie 
h.ive  seen  that  differential  absorption  and  Faraday  rotation  are  related 
to  the  relative  magnitudes  of  the  major  and  minor  axes  of  the  polariza- 
tion ellipse  and  to  the  spatial  angular  rotation  of  the  ellipse  relative 
to  tho  ground-based  antenna.  Consider  a rocket  held  stationary  at  some 
i'oint  .ilong  its  irajertory  such  that  f B and  the  rocket's  tr.inslat  ional 
motion  aro  fixed  but  it  is  still  allowed  to  spin.  Figure  Cl  illustrates 
this  situation.  The  orit-niation  of  the  plane  of  pol  .\r  i zat  ion  (x*^,  v*  ) 
relative  to  the  antenna  plane  (x‘*,  y‘^)  is  determined  by  v’  and  P which 
are  to  be  measured  clockwise  from  the  line  of  intersection  of  the  two 
iilanes.  The  maximum  and  minimum  signals  recorded  over  one  spin  cvcle  of 
the  rocker  are  respectively  proportional  to  the  maximum  and  minimum  pro- 
jections of  K into  the  plane  traced  out  by  the  spinning  antenna.  bet 
Fj  be  the  projei'tion  of  Kj  and  that  of  E^,  where  x*"  is  the 

major  axis  of  the  polarization  ellipse,  is  the  minor  axis,  and 

F.j  .ind  E,  aro  given  by  (14).  Because  of  the  syTTimotry  of  the  problem, 
consider  only  values  of  y'  between  0°  and  180'^.  Let  d^  ,ind  d^  be  tho 
ri'speelive  angular  separations  between  E^'and  E^',  and  the  line  of  inter- 
section ot  the  two  planes.  Since  d^  and  d^  lie  on  opposite  sides  of  the 
line  of  intersection,  d = ^d^j  + itl^l  is  the  angular  separation  between 
F|  and  E,,  . In  gener.al,  the  project  ions  of  two  orthogon.il  vectors  like 

r ,1* 

EjX  and  E^y  of  Figure  C3  onto  a plane  that  is  not  p.irallel  to  the  plane 
ol  the  vectors  (8^0  in  Figure  C'3)  ;ire  not  orthogonal.  If  6 = 71/2, 
then  d = 0 and  jlj  and  E^  are  parallel.  The  angle  between  and  its  pro- 

jection Ej’  is  called  b^. 


K,*  = F,jCos(b, ) 


(C14) 


Similarly  for  E^,, 
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= E2Cos(b2)  (C15) 

Napier's  rules  for  right  spherical  triangles  give  the  relationships 

sin(b|)  = siny'sinB  (C16) 


anil 


sinfb^)  = cosy'sinS  . (C17) 

The  angles  ti  ^ anti  in  the  plane  of  the  antenna  are  then  given  by 

sin(dj)=tan(bj)cot3  (C18) 

sln(d^)  = tan(b^)cotB  • (C19) 

Some  arbitrary  vector  in  the  polarization  plane  E can  be  written  in 
terms  of  the  orthogonal  components  along  and 

E = E + E v*^ 

— X—  y" 

These  components  can  be  projected  onto  the  plane  of  the  antenna 


E = E cos(b,) 

X X 1 

E = E cos(b  ) 

y y 2 

Tile  vector  sum  of  E^^  and  in  the  antenna  plane  is  given  by 

Ie'P  = E'  • E'  = E^cos^b,  + E^cos^b,  + 2E  E cos  (b , ) cos  (b, ) cos  (d ) . (C20) 

''  __  1 y 2 Xy  > '2''''' 


Since  E lies  in  tlie  polarization  plane,  the  magnitudes  of  E and  E are 
~ X y 

less  than  or  equal  to  the  major  and  minor  axes  of  the  ellipse  E and  E 

I 2 

respectively.  Since  d is  not,  in  general,  equal  to  90°,  E^^ and  are 
not  orthogonal  vectors.  Consequently,  may  not  be  the  largest  value 
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pro  j fc  t I’d  into  the  plane  of  the  antenna,  nor  the  smallest.  As  the 
phasor  K rotates  in  the  polarization  plane,  the  minimum  value  of  I E' I 
in  the  antenna  plane  from  (C20)  occurs  when  = 0 . The  components 

E^  and  E^  in  the  plane  of  polarization  satisfy  ^ 


E ^ 

X 


E 2 
1 


(C21) 


Substituting  (C21)  into  (C20) , taking  the  derivative  with  respect  to 
and  setting  the  result  to  zero  gives 


E 

y 


(C22) 


where 


E^  cos(b, ) 


R Ej cos (bj ) cos (d) 


Ej  cos(bj ) 
E20os(b2 )cos(d) 


(C23) 


Rearranging  (C22)  into  the  binomial  form 


gives  the  roots 


1/2+  (1  - R)''^/2 


(C24) 


where 


R 


(C25) 
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*^**‘-’  smaller  value  of  (Ey/l’,)^  corresponds  lo  the  major  axis, 
while  for  Ej  < E^'  the  larger  value  corresponds  to  the  major  axis.  The 
data  becomes  extremely  difficult,  if  not  impossible,  to  reduce  whenever 
E^  < E^  . After  this  situation  occurs,  reduction  of  the  data  will  be 
very  error-prone.  Hence  a discussion  of  this  situation  will  be  neglected. 
Let  Ej  and  E^  be  the  largest  and  smallest  field  components  respectively 
in  the  antenna  plane.  Figure  C3.  By  combining  (C24) , (C21)  and  ((120)  for 
the  usual  case  where  E " ^ E * , 

i 2 

E'"  = (1  + R)  E'V2  + (1  - R)  e;V2 

+ (4  - R^)-  Ej'E^'cosCd)  (C26a) 

and 

Ej^  = (1  - R)  EjV2  + (1  + R)  EjV2 

- (4  - R^)  ^ E|'  E/  cos(d)  (26b) 

The  spin  introduced  error  in  Faraday  rotation  can  now  be  found  by  deter- 
mining the  angular  difference  between  the  "apparent"  minor  axis  E^  and 
the  "true"  minor  axis  pn'Jeccion  E^'  . Since  the  po  1 ar  i.-.at  ion  ellipse 
projected  onto  the  antenna  plant'  is  an  ellipse  also,  the  anp.ular  dif- 
lereni'e  T.  is  given  liy 

= ■siir'|[(Ej/E;)'  - i]  [(e;/e;)'  - iT'l'"  • (e27) 

The  error  in  differential  absorption  is 

r,  = [(e;/e^'-  i)(e;/e;  + d"']  - a (c28) 

where  A is  given  bv  equation  (20)  of  Chapter  II.  For  the  usual  case 
where  ^ F.^  , f,  will  change  sign  as  Y increases  past  a multiple  of 
00  . Similarly,  . will  change  sign  as  Y increases  past  a multiple  of 
somi'  ang,le  between  0°  and  90°. 
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Tlu>  plots  in  Figure  C4  correspond  to  the  case  of  a 10  MHj,  wave  being 
recorded  for  10  seconds  (over  7.4  km)  after  propagating  through  the  summer 
PCA  medium  to  an  altitude  of  77.8  km.  At  this  point,  a differential  ab- 
sorption of  about  -15  dB  has  been  induced  on  the  signal.  The  sole  differ- 
ence in  input  between  the  two  is  that  in  the  upper  plot,  the  electron 

density  has  been  set  to  zero,  while  in  the  lower,  N is  identical  to  the 

e 

PCA  profile  throughout.  Rotational  parameters  of  a = 10  sec  , Y = 0.1  sec 
and  B = 30°  (see  Figure  Cl)  were  taken  as  constants  over  the  second  interval 

in  the  same  fashion  as  indicated  in  (19)  of  Chapter  11,  = S.E  ^ and 

V.,  = £,K^  are  the  recorded  signal  borders  if  there  were  no  coning  present 
(see  Fii'.ure  6B)  , and  V = £[•■  * and  V = tF,  'are  the  distorted  signal  borders 

with  coning  given  by  (C2b)  . Thus  in  the  upper  plot,  the  modulation  of 

and  is  illustrated  by  V^'and  respectively  over  one  complete  cvcle  of 
> . The  slii’dit  decrease  in  the  signal  amplitude  in  Figure  C4  arises  from 
the  increasing,  distance  of  the  rocket  from  the  source.  In  thi'  lower  plot, 
tile  F.iraday  rotation  and  differential  absorption  effects  are  included  to 
illustrate  aitual  errors  that  would  he  incurred  in  the  true  signal  enve- 
lope upon  ignoring  coning,  effects  in  data  reduction.  Note  that  y'  = yt  - F 
no  longer  completes  one  cycle  as  a result  of  Faradav  rotation. 

Tile  coning  errors  induced  into  the  unreduced  data  are  df,/dz  and 
d',/dz.  Thesi'  errors  may  he  estimated  numerically  hv  calculating  vari- 
ations in  (C27)  and  (C'’8)  over  small  distances  in  the  coded  model.  For 
typical  coning  parameters,  when  A is  less  than  -10  dB,  df./dz  can  easily  he 
of  till'  same  m.agnitude  as  dA/dz  or  greater.  The  term  d(i/dz  is  generally 
much  smaller  than  dF/dz  but  can  also  become  significant.  In  general, 
d',/dz  and  d(',/dz  are  functions  of  y',  8 and  Kj/K^.  The  dependence  of  dt;/dz 
and  dy/dz  on  thosi'  variables  is  non-linear  and  aperiodic  with  altitude. 
Siatistii’al  avi'raging  of  the  unreduced  data  may  still  leave  significant 
errc'r . ior  a higher  degree  of  accuracy  it  is  necessary  to  monitor  the 
rocket  oriiTitation  in  flight  continuously  and  develop  a reduction  scheme 
that  will  remove  the  error.  The  remainder  of  this  appendix  will  be  the 
development  of  such  a scheme. 

In  the  plane  of  the  antenna,  tlte  electric  field  E' , is  elliptical 
and  satisfies 


TIME  (SEC) 


TIME  (SEC) 


C-14 


(C29) 


whore  ami  E^  are  the  orthogonal  components  of  e'  along  e'(x’  direction) 
and  E^  (y’  direction)  the  major  and  minor  axes  of  the  ellipse  respectively. 
The  field  components  e'  and  k'  can  then  he  broken  into  the  (x' , ^')  com- 
ponents by 


E'  = E|'cos  (Vj.x'  + Ej'sin  ^y'  (C30a) 

E'  = E/cos  '^x'  + E,'  sin  O^y' 

where  0^  and  0^  are  the  angles  measured  from  x'  to  E^' and  E^' respectively, 
sec  ligiire  C5.  Note  that  d = 0^  + 0^,  where  d is  given  hy  the  sura  of  Idjj 
and  jd^i  in  (CIS)  and  (C14).  Substitution  of  (C30a)  into  (C29)  gives 

= (E'VE_^=^sin^0^  + cos^0j)"‘  (C31a) 

and  (C30b)  into  (C29)  yields 

Ej^/Ej^  = (E3VEJ^sin^02  + cos^CO^,)"'  (C31b) 

Combining  (C26a)  and  (C26b)  gives 


(E^)^  + (eJ)‘  = (E')^  + (Ep^  (C32) 

which  with  (C31)  and  some  rearranging  gives 


Solving  the  above  expression  simultaneously  with  d = |dj  + |d^i  leads  to 
a binomial  form  having  roots 


tanOj=  -[(K  - l)tan(d)]/2  - [[(K  - l)tan(d)]^/4  - K 


(C33a) 
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Figure  C5.  E-fielcl  components  in  tlic  plane  of  the  antenna. 

and 

tan02=  -[(K  - l)tan(d)]/2  - [ [ (K  - l)tan(d)]V4  - K | '^  (C33b) 

where  K = (E'/Ep‘^  . If  d = j d^  | + | d^ ) is  found  from  (C16)  , (C17)  , 

(C18)  and  (C19)  and  K is  known,  0^  and  0^  can  be  found  from  (C33).  Solu- 
tion of  (C30),  (El A)  and  (CIS)  tlu>n  gives  Fj/E^  which  is  the  ratio  needed 

for  the  t Iculation  of  A.  The  rate  of  change  of  A with  altitude  — 

dz 

follows  immediately. 

Analytically,  the  Faraday  rotation  coning  error  introduced  by 
measuring  the  difference  between  signal  minima  and  spin  marks  is 

d z A z d z 


L. 


(C34) 
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where  AF'  is  the  ohan^te  in  anj;iilar  separation  between  signal  Tnlnima 
and  spin  marks  over  the  interval  Az,  and  dF/dz  is  the  true  Faraday 
rotation  rate.  The  calculation  of  dlj/dz  is  not  difficult,  and  then 
dl’/dz  can  be  found  from  the  measured  value  AF'/Az. 

Consider  an  altitude  interval  Az^  t'lat  is  small  enough  that  the 
(|uantities  dl'/dz.  and  dA/dz.  may  he  considered  constant  over  this  inter- 

■J 

val . These  quantities  uniquely  determine  electron  density  and  colli- 
sion frequency  within  this  interval.  This  will  permit  reduction  of 
the  data  measured  over  Az.  independent  of  measurements  at  other  altitudes. 
Ihe  ratio  of  the  minimum  to  [ruiximura  signal  as  shown  in  Figure  9 is  the 
quantity  ~ ^issuming  linear  reception.  The  parameter  F^ 

is  thi'  angular  separation  between  signal  minima  and  spin  marks.  The 
angular  separat  itni  between  precession  marks  and  signal  minima  is 
'i'.  - yt  - Y..  There  is  an  error  term  in  this  assumption,  arising  from 
the  presence  of  i.  in  F . , as  in  (C34) . This  is  ultiirnttely  a second- 
order  error  term,  however,  since  y . is  used  for  correcting  coning  error 
itself.  The  precession  marks  can  be  determined  from  magnetometer  data 
by  a process  similar  to  determination  of  the  spin  marks.  Angle  6^ 

I an  be  determined  from  the  amplitude  of  the  magnetometer  data  at  z.. 

7. 

With  the  lour  tabulated  lunctions  .,  y , F^,  and  p.  discussed 

above,  the  data  rediu-tion  for  the  k^*'  interval,  Az,  = z,  - z , is: 

k k k-i 

n 

2) 

3) 

4) 

5) 

h) 

7) 


C.ilculate  from  (Clb),  (Cl  7)  , (C18),  and  (C19)  the  parameters 
P',  . b^ , d;  , d'  , and  d.  = d;  + d^  for  i = k,  k-1 . 

Calculate  G'^  and  0^,  f = k,  k-1  from  (C31). 

From  (C31)  the  ratios  (Ej/Ep^|  . and  (E'/Ep^|^.  , i = k,  k -1 
can  be  calculated. 

Calculate  K /E  I . from  (C14)  and  (C15)  for  = k,  k-1. 

12': 

From  (18)  and  (17)  calculate  A:  for  i = k,  k-1.  Then 

dA/dz;  = (A.  - A.  )/(z , - z,  ■'  ) 

k k k-1  k k-1 

Calculate  f,,-  from  (C27)  for  i = k,  k-1.  Note  that 
= [^'/E;|.][E;/E;i..]->  . 

Calculate  dF/dz],  from  (C34)  where  AF,'  = F'  - F'  , 
k k k k-i 
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a) 


Iterate  collision  frequency  V and  electron  density  N bv 
, . ■ m e 

the  expressions 


dF/dz||^. 

, dF/dzL  ^ dA/dzl,  v'^ 

^ .__k  k m 

m dlVd^l^  dA/dz!|^. 

k " k^  I I 

where  N and  V are  initial  assumed  values,  and  dF/dz  , ^ and  dA/dz  , ^ 
e m , k k 

K K 

are  the  rate  e<iuatlons  found  by  subst  i tut  inc,  V and  N into  the 

me 

Sen-Wyller  equations.  By  replacing  k*  with  k,  the  above  expressions 

k k 

rapidly  converge  to  N and  Y • This  last  step  is  used  bv  r1Pc‘kt7,u  ct  al . , 

e m ■ ^ 

[1967]  and  Bennett  et  at.,  [l972]. 

It  was  shown  in  Chapti-r  III  that  reliable  differential  absorption 
d.ita  can  be  obtained  at  somewhat  lower  altitudes  than  Faraday  rotation 
data.  At  these  low  altitudes,  coning  errors  are  insignificant.  The 
procedure  given  by  Mcchtlii  ct  al . , [l967]  and  Bennett  ct  al.,  [l972]  is 
to  extrapol.ite  the  collision  frequency  to  lower  altitude  by  assuming  it 
is  proportional  to  pressure;  = k p(z).  Piakenson  et  al . , [l976] 

state  that  a propot ional ity  to  density  P(z)  may  be  a better  assumption, 
however.  The  first  equation  in  step  (8)  above  can  be  used  to  determine 
the  lower  electron  density  values  if  is  extrapolated. 

In  actual  practice  there  may  prove  to  be  better  ways  to  determine 
the  quantities  F^,  Y , and  6 but  this  is  beyond  the  scope  of  this  paper 
and  will  not  be  discussed. 


